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Corrections to listing of CANDO: ERL-M251 


SUBROUTINE PT: 
Card C-246 should read 
IF (KKTT, EO.1) go to 62 
Card C-253 enone. be replaced by the following three cards: 
I= wm-l 
DO 35 J = NN, NB 
= 
The following card should be added right after C-257 


62 it = i aoa 


py a2) 


itl 


ERL-M251 ERRATUM: 


i+ 
Page 37, the expression for Hq i should read: 





. - y : : : : see : 1 
s* * (s*) txt * (g'**-g')] * [H * (g =e") 
a pee ET, Sa 
(Gg) 2 (ES) (gttg') * a * (g =e) 


Page 38, the indicated blocks should contain the following 





Compute the following scalars 






(ey See) 


ey, saad mus 
(g -g) *W *(g -g) 


~ 







Compute matrices 











Saree rs eer 
[x *(g -g)] * [H * (ge )] 


~ 
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ABSTRACT: 

This report describes the operation and use of CANDO (Computer Analysis 
of Networks with Design Orientation), a time domain, design oriented, 
analysis program for linear time-invariant networks. The networks may con- 


tain dependent and independent sources, capacitances, resistances and 


‘inductances. 


The network analysis problem is to obtain the branch currents and 
voltages by solving a set of simultaneous differentio-algebraic equations 
derived from the branch relations and Kirchhoff's current and voltage laws. 

In the program CANDO, nonindependent source tree voltages and non- 
independent source link currents form a basis set of variables. This formu- 


lation yields the automatic satisfaction of Kirchhoff's laws and also allows 


us to solve an optimally low order system of equations. A numerical integration 


formula reduces the system to a set of simultaneous algebraic equations of 


the form 


which are solved by minimizing a performance criterion € defined by 
a 2 
C= 2 ) 8 
i | 

CANDO utilizes the Fletcher-Powell and Rohrer search minimization 
algorithms in minimizing ¢. 

The tree-picking and internal current scaling algorithms are such that 
large value spread and large time constant spread problems can be handled 
reasonably effectively and efficiently. 


The program CANDO is written in FORTRAN IV for the CDC 6400 system. 
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INTRODUCTION: 

In computerized fully automated network design, where one performs 
seyerel network analyses at each iterative step of the design process, 
one has a need for an analysis program which can be made arbitrarily 
accurate, and whose execution time is a function of the accuracy desired. 

CANDO satisfies the above requirement, since one can specify either a 
maximum permissible value of the performance criterion ¢, and/or a maximum 
number of minimization iterative steps, per time point. 

The iterative scheme utilized by CANDO is such that it can be easily 
extended to networks containing nonlinear and time-varying elements » by the 
simple addition of corresponding subroutines. 

Section I describes the iterative scheme utilized in solving both the 
initial condition problem and the general time point problem. Section II 
describes the MAIN program and subroutines associated with CANDO. The 
appendices describe the use of CANDO in solving network problems, and the 


way in which dependent sources must be modeled. 
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SECTION I: Network Analysis 
Theory : 


The solution of a network consists of finding the vectors Lb and Xy 
representing all of the branch currents and voltages, respectively. The 


Solution is obtained by solving the following simultaneous set of equations 


£0, ; i,) = 0 branch relations 


con =O Kirchhoff's current law 
By, = 9 Kirchhoff's voltage law 


where Q is the fundamental cutset matrix and B is the fundamental loop matrix, 
Upon the selection of an appropriate tree, we may renumber the NB branches of 
our network in such a way that the lst NN-1 (where NN is the number of nodes 
in our network) branches form the tree. 


With this numbering scheme » We may partition Q and B as follows: 
Q=(ZiF] and B= [-F'} 7] 


where I is the identity matrix and hence we have, from Kirchhoff's laws, 


4, => BH, 
= fF 
yy =F > 


where the subscripts t and # refer to tree branches and links respectively. 
Our tree-picking algorithms are discussed in Subroutine PT; suffice 

it to say that in both our tree-picking schemes, independent voltage 

sources must be tree branches and independent current sources must be links. 


Upon using a numerical integration formula (see Subroutine ERROR), our set 


'® 
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of linear equations reduces to 


E(v,» i,) = 0 
and the dimension of the system is NB minus the number of independent 
sources. 


This set of equations is solved by minimizing a performance criterion 


€ defined by 


i 21, 
e-3) & =3E-E 


The minimization is performed by using the Fletcher Powell algorithm 
to give us a direction along which to minimize (see Subroutine FLPOW), 
and using Rohrer search to find the optimum step size along that direction 
(see Subroutine FNDALFH). 

This iterative scheme is continued until either the error is suffi- 
ciently small or we have gone the desired number of steps. (See Program 


MAIN). 


Solving the Initial Condition Problem: 

Given the initial capacitance voltages and inductance currents, we 
wish to compute My and i, the branch voltages and currents at the initial 
time point. This reduces to the problem of solving a coupled system of 
algebraic equations. 

A proper tree is selected according to the algorithm described in sub- 
routine PI, allowing us to reduce the order of the system. This algorithm 
maximizes the number of capacitances in tree branches and the number of 
inductances in links. In general, if there are no capacitance or capaci- 


tance-voltage source loops, or inductance or inductance-current source 


(a 


cutsets, all capacitances will be tree branches, and all inductances will 
be links. 

Kirchhoff's current and voltage laws become automatically satisfied 
upon the selection of the tree voltages and link currents to be a basis set 


of variables, with the following relations emerging: 


x= 7 Ey 
and 
‘at t 
My =FM > 


where the subscripts £ and t refer to links and tree branches respectively, 
and where Q = [I ,F] is the fundamental cutset matrix based on the proper 
tree. 

To further reduce the dimension of the system, we do not treat independent 
sources as variables, but rather require independent voltage sources to be 
tree branches, and independent current sources to be links. 

Furthermore, tree capacitances are treated as independent voltage 
sources with their voltage being set to the initial condition. Similarly, 
link inductances are treated as independent current sources with their 
current being set equal to the initial condition. Tree inductance voltages 
and link capacitance currents are set equal to zero, and their initial 
conditions are lost. Hence, incompatible initial conditions (e.g., in- 
compatible capacitance voltages around a capacitance-voltage source loop or in- 
compatible inductor currents at an inductor-current source cutset) are forced 
to be compatible, with the tree inductances and link capacitances loosing 


their specified initial conditions. 
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Thus the system of equations we are solving is one corresponding to 
a network consisting only of independent sources, dependent sources and 
resistances. 

The solution is obtained in a way analogous to the solution of the 
equations associated with the "general time point problem" (see next 


section), and will not be detailed here. 


Solving the general time-point problem: 

Knowing i, and Yy> the branch currents and voltages at the previous 
time-point, we wish to compute i, and Xos the branch currents and voltages 
at the succeeding time point. The problem reduces to solving the coupled 
system of differential and algebraic equations, consisting of the branch 
relations, and Kirchhoff's current and voltage laws. 

A tree is selected according to the optimal tree algorithm (see sub- 
routine PT), allowing us to reduce the order of the system, and hence to 
improve the convergence of our iterative minimization scheme. 

The network is renumbered such that the first NN-1 branches form the 
tree, and the subsequent branches are links. Kirchhoff's current and 
voltage laws become automatically satisfied upon the selection of the tree 
branch voltages and link currents to be a basis set of variables. 


The tree currents i. 4 are given, in terms of the link currents i, by 


i,°=- B, ’ 


where Q = [Ii F] is the fundamental cutset matrix associated with our 


optimal tree. 


Similarly, 


where F is the transpose of F. We thus wish to solve the set of equations 


? B(i,, x) = 2 
2 which are our branch errors (see subroutine ERROR). 

To further reduce the dimension of the space, we do not treat independent 
source values as variables. Notice that this requirement dictates that inde- 
pendent voltage sources be tree branches, and that independent current 
sources be links. 

The solution to E(i,, x) = i) is obtained by minimizing a performance 
criterion ¢€, given by 

€ =) BE,” 
i 

Note that we have an exact solution when « = 0. 

€ is minimized via the Fletcher-Powell minimization algorithm‘) 

(see subroutine FLPOW) with Rohrer search (see subroutine FNDALPH) being 
utilized to find the optimal step size along the direction generated by 
the Fletcher-Powell algorithm. 

Once % and i, are known, we have 

-~ F 
dg *|-- 7 
x 
£ and 
i 
. Xe * | 7 xe? 


and hence, we proceed to the subsequent time point in exactly the same way. 
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SECTION II: PROGRAM C 


MAIN program: 

The MAIN section of program CANDO ensures that the desired subroutines 
are called in the correct sequence. It solves both the initial condition 
problem and (when required) the general time point problem, in the way 
discussed in Section I under Theory. It also controls the calls to the 


output subroutines. 


os 


BLOCK DIAGRAM: 


MAIN PROGRAM 


CALL READIN 


- Read in network description 


- Compute proper tree (call PT) 


- Compute associated F natrix 
(call FCSM) 


CALL KONST 
- Compute branch element constants 


- Select internal current scale 
(if desired) 


- Initialize tree voltages and link 
currents at initial time point 


CALL INCREM (May call PERIOD) 


Initialize nonconstant sources 
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CALL CALCAL 


Compute tree currents and link 
voltages from tree voltages 

and link currents via Kirchhoff's 
laws 





Have 
we exceeded 
desired mmber 
of iterations 






CALL desired 
output 
routines 


CALL ERROR 














- Compute branch errors 


- Compute performance 
criterion 










Is 
performance 
criterion less than 
specified 
error 


CALL. FNDGRAD, Compute. gradients 


CALL FLPOW, Compute minimization 
direction 


CALL FNDALPH, Compute optimal step size 


CALL CNGVARS, Change variables 


‘we 










Have Go to next 
we gone far problem to 
enough in time be solved. 





CALL INCREM (which may call PERIOD) 






- Compute independent source values 
at next time point 











Is 
a new 
tree required 





- compute optimal tree 






CALL FCSM - compute associated 
F matrix 






CALL KONST - compute new internal 


scale (if desired) 






compute new branch 
constants. 
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SUBROUTINE READIN: 


This subroutine reads in the complete network description along with 


all pertinent output and control specifications. As a data card check, 


- the given network description and the pertinent control specifications are 


printed out. 

For a detailed description of the variables read in, see the section 
entitled "Data Cards for CANDO." 

This subroutine also calls subroutine, PT which computes an appropriate 
tree to solve the initial condition problem. Te no independent voltage 
source loops or independent current source cutsets exist, subroutine FCSM 


is called, otherwise the program proceeds to the next set of data. Sub- 


routine FCSM computes the nontrivial portion of the fundamental cutset 


matrix. 


.s 


BLOCK DIAGRAM: 


Input and output network description 
and associated pertinent data. 


CALL PT. KLOOP will be set to 1 if 
there is either an independent 
voltage source loop or an inde- 
pendent current source cutset. If 

KLOOP = 1, the program proceeds to 

the next network to be analyzed 

















SUBROUTINE READIN 
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SUBROUTINE PL 

The first time it is called, this subroutine computes a proper tree 
for the initial condition problem; on subsequent calls in the same problem, it 
computes an optimal tree. The distinguishing features of the above two 
tree-picking schemes will be enunciated below. 

If the given network has NN nodes and NB branches, it will have NN-1 
tree branches and NB - NN + 1 cotree links. 

A proper tree consists of the first NN-1 branches which do not form 
loops, with the following selection priority scheme (in descending order ) 

1. independent voltage sources 

2. dependent voltage sources 

3. capacitances 

hk, resistances in ascending order of value 

5. inductances 

6. dependent current sources 

7. independent current sources 

An independent voltage source appearing in a link indicates the 
presence of a loop of independent voltage sources. An independent current 
source appearing in a tree branch indicates an independent current source 
cutset. In both of these cases, an error diagnostic will be printed out, 
and the program will proceed to the next problem. 

An optimal tree consists of the first NN-1 branches which do not form 
loops, with the following selection priority scheme (in descending order) 

1. independent voltage sources 

2. equivalent resistances in ascending order of value 


3. independent current sources 


‘g 


a 


The equivalent resistance values are given by the following prescription: 


Controlled voltage source ~ 107° 


Resistance (R) > R 
Capacitance (C) > at 
Inductance (L) eo oe 

2H 


' Controlled current source ~ 10°, 


where H is the specified integration step size for use with the trapezoidal 
rule. If C = 0, the corresponding equivalent resistance is set to 10°", 

The internal branch numbering scheme is such that the first NN-1 
branches are tree branches, the remaining being cotree links. This modification 
allows us to subsequently store only a portion of the fundamental cutset 
matrix. 

Note that in the selection of a proper tree, the independent current sources 
are placed at the end of the branch list. In the subsequent selection of 
an optimal tree, the independent current sources are not moved around and 
hence need not be rearranged. 

This subroutine, when desired, prints out the selected tree and the 
rearranged network. 


Note that when a network contains only sources and resistances (i.e., 


no reactive elements), the proper and optimal trees coincide. 


vd 
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BLOCK DIAGRAM 
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Reorder branches according 
to desired priority scheme 
(proper tree or optimal tree) 


















Select first NN-1 branches 
which do not form loops to 
make up desired tree 








Are there 
independent 
voltage source loops 
or independent 
current source 







Print error 
diagnostic 


Read in next 
set of data 












If required, rearrange link 
branches placing independent 
current sources at the end 

of branch list 


SUBROUTINE PT 
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SUBROUTINE FCSM (Fundamental Cutset Matrix): 

Given a connected network, and a specified tree, there exists a unique 
(NN-1 x NB) fundamental cutset matrix Q, where NN and NB are the number of 
nodes and the number of branches in the network, respectively. 


Kirchhoff's current law is then given by: 


5 


where i, is the branch current vector. 
If the branch numbering is selected such that the first NN-1 branches 
are tree branches (see subroutine PI), then, Q can be partitioned as 


follows: 
Q = (ZI) 


and only the (NN-1 x NB-NN + 1) F matrix need be stored, resulting in a 
substantial saving of computer memory. 

With each link, of our connected network, is associated a unique funda- 
mental loop, consisting of the link itself and sufficient tree branches 
required to close the loop. Each column of F corresponds to such a funda- 
mental loop. This intermediate scheme is used to compute F. 

If desired, this subroutine can output the F portion of the fundamental 


cutset matrix. 


BLOCK DIAGRAM: 


Select a link 







Search along tree branches 
until the corresponding 
unigue fundamental loop 

is obtained. 


Fill the corresponding 
column of F 





SUBROUTINE FCSM 
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SUBROUTINE KONST: 

This subroutine computes the branch element constants required in the 
branch relations. It initializes tree voltages and link currents, and, when 
desired (NSCALE # 1), it selects a new current scale factor each time a 
new tree is picked. 

There are no constants associated with independent sources. | 

The constants are computed according to the prescription given below, 
where H denotes the integration time step, in seconds, 8S denotes the scale 
factor which selects the desired current unit (e.g., S = 1000 implies current 
unit is milliamps), and A, denotes the branch element constant associated 
with the ith branch. See subroutine ERROR for the form of the branch rela- 
tions. 


Dependent voltage sources: 


a) voltage controlled case 


A, oK, 


where K, is the ith coupling constant and is unitless. 


b) current controlled case 


~ 
n 
a 


where K, is the ith coupling constant, in ohms. 


Capacitances: 


a) tree branch case 


iv 


b) link case 





* ¥ 
— 2 C5 S 
i “HH | 
where Cc, is the ith capacitance value, in farads. 
Resistances: 





a) tree branch case 


Ry 
Fae 
b) Link case 
=o 
Aa 
1 


where R, is the ith resistance value in ohms. 


Inductances : 
a) tree branch case 


2*L 





et ee 
A, eS 
b) link case 
H*S 
“BFL, 


where L; is the ith inductance value in henrys. 


Dependent current sources: 


a) voltage controlled case 


= x* 
A, K, iS) 


where K, is the ith coupling constant, in ohms, 


b) current controlled case 


where K, is the ith coupling constant and is unitless. 


Initialization; 

The first time that subroutine KONST is called it initializes tree 
voltages and link currents as follows: 

tree voltages: Tree inductance and tree dependent current source 
voltages are set equal to zero, the remaining being set to the given initial 
conditions. 


dink currents: link capacitance and link dependent voltage source 








currents are set equal to zero, the remaining being set to the given initial 
conditions. 
It should be noted that the above initialization scheme makes sense 


in light of the fact that the associated tree is a proper tree. 


Automatic computation of scale factor: 
When desired, this subroutine will compute a new current scale factor 
each time a new tree is picked. 
The algorithm, for selecting the optimum scale factor, for both proper 
and optimal trees, is as follows: 
1. select smallest link resistance value (say R,) 
2. set internal current unit to be a. where R, is in ohms, e.g., 
if R, = 14600 ohms, the internal current unit will be TT 
milliamps. 
If our tree selection yields the fact that there are no link resistances, 
the scale factor is set equal to the maximm tree resistance. If our network 


contains no resistances, the scale factor is automatically set equal to 1.0. 


SUBROUTINE CALCAL: 
This subroutine computes the tree branch currents i, and the link 
voltagesyv , from the link currents i 2 and the tree voltages My? respectively, 


via 


f 

i] 
| 
pe 


and 


where F is the nontrivial part of the fundamental cutset matrix (see Sub- 
routine FCSM), and where the prime denotes the transpose. 

Thus we see that Kirchhoff's voltage and current laws are automatically 
satisfied, for any given network, because i, and v, are forced to satisfy the 


2 
above relations. 
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SUBROUTINE INCREM: 

Upon the completion of computations at one time point, this sub- 
routine initializes the program for the next time point. 

The most important aspect of the above mentioned initialization, is 
the computation of the independent poets values at the next time point. 

When periodic independent sources exist, INCREM calls subroutine PERIOD 
which then computes the corresponding signal values. 

An arbitrary independent source is specified by two arrays (up to 100 
points each) denoting the time points and the corresponding source values. 
If the program time does not correspond to a time point, linear interpolation 


between the previous and the next time points is automatically carried out. 










e.g. 

independent x--Data Read in for arbitrary waveform 
SQUECE __. Waveform assumed by program 
value 


value assumed 
by program 


Program 
time point 


There are three basic ways in which a nonconstant independent source 


can be specified: 


ee 
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1. arbitrary (as above) 

2. exponential 

3. sinusoidal 
To see how these are read in, see the section on "Data Cards for CANDO." 
The values of the independent sources, beyond the last specified time 
point, are taken to be zero. 
For sinusoidal and exponential source specifications, the program uses 
library SINE and EXPONENTIAL routines. 


The general operation of this subroutine may be seen on the following page. 


ad 


BLOCK DIAGRAM 


Do the following for each nonconstant 
independent source 













Call subroutine 
PERIOD | 






exponential 
type 


Compute source 
value explicitly 


sinusoidal 
type 
















fompute source 
value explicitly 










Interpolate, between 
adjacent time points, 
linearly, to obtain 

source values 







SUBROUTINE INCREM 
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SUBROUTINE PERIOD: 
This subroutine computes the signals of periodic independent sources, 
for the next time point. 


Linear interpolation is used, as in subroutine INCREM, the only dif- 





ference being that when the last specified time point is exceeded, the 


waveform is repeated. 


e.g. 


waveform read in 


t 





waveform generated by computer 


Program 
time points 


Thus we note that when the last specified time point does not corre- 
spond to a program time point (i.e., t) + 0H where n = 0, 1, .--), the 


period, T, of the waveform is increased so that T = nH, n=1, 2, .... 


‘a 
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SUBROUTINE ERROR: 
This subroutine computes the branch errors E, and the performance 


criterion which is to be minimized, defined by 


where i ranges over all. nonindependent source branches. 
Note that when both the branch relations and Kirchhoff's laws are 


simultaneously satisfied, ¢€ will be identically zero. 


Branch errors: 

Notation: The A A denote the branch element constants defined in subroutine 
KONST. Vio 
Vaz and C.4 denote the branch voltage and current, at the previous time 


point, respectively. The Ey of course, denote the branch error associated 


with the ith branch. 


integration: 
The trapezoidal integration scheme is utilized to perform the inte- 


grations associated with the reactive element branch relations. 


The integral of a function f over an interval H is approximated by 


t wl 
ie f db =3 Ee + ze 


Since an iterative minimization scheme is utilized at each time point, 
the. trapezoidal rule, in this context, becomes equivalent to a multistep, 
second-order predictor corrector scheme. 

We always integrate either voltage or current. Since both are treated 


alike, only the current integration will be considered below. 


and C 42 denote the present branch voltage and current respectively. 


We have the following typical error term 
=e - * a 
E. (Vio V53) +A, (C3, + C54) 
where Vax and Cia are known. 


As the error is reduced to zero, we have, in effect, 


Nao aa (Cy, + C54) 


i.e., the left hand side is the value of the integral, and Cio is chosen 


so as to satisfy the above equality via linear interpolation. 





Note that the area under the linear interpolation is the same as the area 
under the true curve. 

The above argument indicates that for this integration scheme to be 
accurate, the step size must be much smaller (of the order of 1/10 is 


sufficient) than the smallest contributing time constant. 


Dependent voltage sources 


a) voltage controlled case 


B= - Vag tay * Vo 
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where the subscript K denotes the controlling branch. 


b) current controlled case 


E, =- V3 


x 
i ‘if xe Ke 


2 


where K again denotes the controlling branch. 


Capacitances: 


a) ‘tree branch case 


B= - (Vin - Vay) + Ag * (Cyo * C44) 


b) link case 


BL =~ (Cio +C,,) +a, * (Vi - Vyy) 


Note that at the first time point, when we are solving the initial condition 
problen, E, is set to zero. 
Resistances: 


a) tree branch case 


ty 
i} 


= * 
V; + A; Cc, 


i 2 2 
b) link case 
= « * 
By =~ Gin tA, * Vip 
Inductances: 


a) tree branch case 


B= > (Wyo * V5) + Ag * (Cy - Cy) 


b) link case 


EB, = - (Cyy - Cy) + Az * (Vy + Vay) 


Note that E, is set equal to zero at the first time point, while solving 
the initial condition problem. 
Dependent current sources: 


a) voltage controlled case 


EB, = - Cio tay * Vip 


where K denotes the controlling branch. 


b) current controlled case 


EB, =-C.,+A 


* 
377 Cin TAG * Cy 
where K denotes the controlling branch. 
Performance criterion: 
The performance criterion € is given by 
€ =) 2? 
i 
i 


where i ranges over all the nonindependent source branches. 
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and from (1) we have, 


Be 
a 


OE 

oy, 
mf 
2 oe 

oy, 


OE 
ge = (rF] ON 
av, I SLE ta 

Of 

oy, 

which reduces to 
at pr ab 
Link current gradient vector: 
Since 
oe 
abo i, 
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(3) 


(4) 


(5) 


(6) 
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we may proceed exactly as in the tree voltage gradient case, and arrive at 


de of OE 
aS -F'- (7) 
dip Oi, * © Oh, 
Branch Gradients: 
OE OE o£ 
The branch gradients ss and 3 7— are computed as follows, where a 


ana S- are branch gradients associated with the ith branch: 
i 


Tree branches: 
Independent voltage sources: 


a) not associated with controlled sources 


OE OE 


Dependent voltage sources: 


35 


Dependent current sources: 


wy, =0 
Oe EB 
5 i 
Remaining tree branches: 
i =), * 5, 


Link branches: 
PPP PPD PPP 
Independent current sources: 


a) not associated with controlled sources 


b) associated with controlled source K 


oe _ 

ov, = Ay * By 
oe _ 

oi. 0 


L 


Dependent current sources: 


Dependent voltage sources: 


Remaining link branches: 


Note, when solving the initial condition problem, the gradients with 
respect to reactive element currents and/or voltages, are set equal to 


zero. 
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SUBROUTINE FLPOW: 

This subroutine computes the direction § along which we wish to minimize 
the performance criterion e€. 

The Fletcher-Powel1‘) aigoritim is utilized with the initial direction 


being that of steepest descent (i.e., along the gradient). 


Fletcher-Powell Algorithm: 
Let H be a positive definite symmetric matrix. 


The direction along which we perform the minimization is given by 
where the superscript i indicates that we are at the ith iterative step, 
and g is the gradient vector. 

~ 


H is updated as follows 


H=I 


itl ; m s* * (s*) Ce i s*)H(g "+ s*)' 
~  *ghyt ® (gt ~ st) (gi - sty eg - 8°) 


Thus we see that use of the Fletcher-Powell algorithm requires us to store 


ann Xn matrix, the Hessian matrix. 


Block Diagram: 





Is 
this the 
first iteration 
for this time 


Compute the following 
scalars 


(s*)' * (eo a s*) 


aie . yi? ‘ s*) 









s* 


~ 














Compute matrices 
: ' 

itl i itl i 

g - SEG" - 8) 

Compute Ee 


act git? = pitt itl 












SUBROUTINE FLPOW 
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SUBROUTINE FNDALPH 
This subroutine computes a » the optimal step size at the ith iterative 
step, which minimizes ¢ along s", a direction computed by subroutine FLPOW. 
The constant ot is such that 
min e(x™ - As”) = e(x" - a's") 
A 
where x is the set of nonindependent source tree voltages and link currents, 


at the ith iterative step. 


The analytical expression for a* is: 


(g*)' * st 


e(x - s+) + (g)' # 8° - e(x) 


i 
a = 


mle 


where g is the gradient at the ith iterative step. 


The situation may be depicted graphically as follows: 





SUBROUTINE CNGVARS: 
Given the direction along which we minimized, st >» and the optimal 
step size ot » the nonindependent source tree voltages and link currents 


are modified as follows: 


ee S je - at * gi 
x2 i, 
new old 


It should be emphasized that only the basis set of variables, i.e., 
tree voltages and link currents, are changed. Tree currents and link 
voltages follow directly from the implicit formulation of Kirchhoff's 


Laws. 
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SUBROUTINE READOUT: 
This subroutine outputs all the desired currents and/or voltages, in 
the desired order, and in the original network numbering scheme. Hence 
for any branche one ‘ae output its current and/or its voltage. 
Note that a one to one correspondence, between desired output variables and 
their internal ordering, needs to be calculated each time a new tree is 


selected. This is done in subroutine PT. 


ke 


SUBROUTINE GRAPH: (Graphical outputs) 

The purpose of this subroutine is to provide graphical outputs of up 
to five variables. For each time point at which an output is desired, the 
main program calls this subroutine. The desired variables are then stored 
in an array. When either the final desired output time point or the two 
hundredth output time point is reached, the variables will be scaled, 
individually, such that maximum use will be made of the output page width. 
The corresponding variable value is also outputed. 


The following figure illustrates a typical graph output, 


Corresponding 
Time Variable Value Branch 8 Voltage 


PA EEE b 
EEPEEEEEEEEL EEE 





MEPEEEE ETT TT EEE EEE 


equal to maximum absolute 
value of signal 


The operation of this subroutine may be seen on the following page. 


For each (up to 5) of desired 
variables, store value 


Is 
the final 
desired or the 
200th output time 
point reached? 












No 


Yes 


Find the maximum absolute value 
corresponding to each variable 
Scale variables according to 
above 


Plot appropriate graphs 


SUBROUTINE GRAPH 
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SUBROUTINE ALLOUT: 

This subroutine outputs all the tree voltages and link currents 
(including those associated with independent sources), in the optimal or 
proper tree numbering scheme. Hence, when it is used, one should also 
output the corresponding tree information which will immediately yield the 


isomorphism between the original network and the new topologies. 
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APPENDIX A: 
Dependent source modeling: 

Each dependent source requires two branches for its complete specifi- 
cation, both of which must be a pare of the network description, and hence 


must be read in as data. 


Voltage controlled voltage source: 


Ideal model: 


oO 
+ 
Vv + V= kv 
Oo———O 
Controlled Controlling 
voltage voltage 


CANDO model: 


Controlled 
voltage 


Current controlled voltage source: 


Ideal model: 


Controlled 
voltage 





CANDO model: 


Vv : v=o 


Controlled Controlling 
voltage current 


Voltage controlled current source: 


Ideal model: 
oO 
+ 
I v 
c 
o———o 
Controlled Controlling 
current voltage 
CANDO model: 





e, 


Controlled Controlling 
current voltage 


Current controlled current source: 


I : I 
c 


Controlled Controlling 
current current 





Ideal model: 
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CANDO model: 


ff 
5 





Controlled Controlling 
current current 
Thus we see that in all cases the controlling current is taken to be 
the current through a zero-valued voltage source, and the controlling 
voltage is taken to be the voltage across a zero-valued current source. 
The coupling constants, denoted by k, are those that should be read 


in as corresponding elements of the VALUE array. 


APPENDIX B: 

Data Cards for CANDO 

Notation: 
I => Integer format 
E => Exponential or floating point format 
A => Alphanumeric format 


Col = Column on data card 


Card #l: 
Variables read in, in sequential order 
NN, number of nodes (I) 
NB, number of branches (I) 
TSTART, starting time (E) 
TEND, end time (BE) 
H, time increment (E) 
NSTEP, number of time iterations per output (I) 
EPS, error criterion (E) 
NCONT, Tree and F output control (I) 


NSTEP 


1 = outputs desired 


NSTEP = 0 = outputs not desired 


col 1-5, 


NN 
col 6-10, - NB 

col 11-25, - TSTART 

col 26-40, - TEND DATA CARD #1 
col 41-55, - H 

col 56-60, - NSTEP 

col 61-75, - EPS 


eol 80, - NCONT 
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Card #2: 
Variables read in, in sequential order 
NGRAPH, number of graphical outputs (I) 
NALLOUT, control variable for use of ALLOUT subroutine (I) 


NALLOUT = 1 = use of ALLOUT is desired 


NALLOUT = 0 = use of ALLOUT is not desired 
JOUT, number of outputs desired (I) 
To be used in conjunction with subroutine READOUT. 
SCALE, scale factor (E) 
€.g., scale factor of 10° sets current unit to milliamps 
NITT, number of iterative minimization steps, per time point (I) 
(NITT = O if using only error criterion) 
NSCALE, control variable for use of internal, automatic current scaling. 
._ NSCALE = 1 = use scale factor read in as SCALE 


NSCALE #1 => use internal, automatic current scaling algorithm. 


col 1-5, - NGRAPH 
col 6-10, - NALLOUT 
col 11-15, - JOUT 
DATA CARD #2 
col 16-30, - SCALE 


col 31-35, - NITT 


col ho, NSCALE 
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Network description data cards 
For each network branch, the following card (or set of cards) is 
needed. The order in which network branches are read in is arbitrary. 
Variables read in, in sequential order 
TYPE, Branch type (A) 
E - independent voltage source 
- controlled voltage source 


- capacitance 


V 

Cc 

L - inductance 
R - resistance 

I - controlled current source 

J - independent current source 
IBRAN, Branch number (I) 

SORTYPE, Independent source type (A) 


C - constant source 


E 


exponential source 
P - periodic source 


S sinusoidal source 


T - time-varying 
CONTYPE, dependent source controlling type (A) 
V - voltage controlled 
I - current controlled 
KONBRAN, dependent source controlling branch (I) 
LEAV, node which branch leaves (I) 
LENT, node which branch enters (I) 
NCARDS, a flag signaling that more data pertaining to this branch needs 


to be read in (independent sources only). (I) 
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NCARDS = 0 = no more data needed 


NCARDS # O => more data needed 


VALUE, value of branch element (E) 
Resistances in ohms 
Inductances in henrys 
Capacitances in farads 
Dependent source coupling constants in ohms, mhos, or unitless. 
COND, - for inductances, the initial current 
- for capacitances, the initial voltage 
- for constant independent sources, the source value 


(in volts or amperes). (5) 


col l, - TYPE 

col 2-4, - IBRAN 

col 5, - SORTYPE | 

col. 6, - CONTYPE DATA CARD FOLLOWED BY OTHER 
col 7-9, - KONBRAN RELATED CARDS IF NCARDS # 0 
col 11-12,- LEAV 

col 14-15,- LENT 

col 17-18,- NCARDS 

col 21-35,- VALUE 


col 36-50,- COND 


If NCARDS # 0, we need the following card(s), to describe the corre- 


sponding nonconstant independent source. 


- If SORTYPE = S, i.e., a sinusoidal source, the signal is assumed 
to be of the form 


A sin (w * t - 0) 


one card is required to describe the above 
col 1-10, signal amplitude A, in volts or amperes (5) 
col 11-20, signal frequency w, in radians/sec (E) 
* col 21-30, signal phase >, at initial time (TSTART), (E) 
~ If SORTYPE = BE, i.e., an exponential source, the signal is assumed 


to be of the form 
aev(*-%) 


one card is required to describe the above 

col 1-10, signal amplitude A, in volts or amperes (£) 

col 11-20, time constant y (sec’+), (E) 

col 21-30, signal phase >, at t = TSTART, (E) 

- If SORTYPE = T or P, i.e., a time varying or periodic source, we 
need the following set of cards, 

a) card #1, col 1-5, number of time points (I) 

b) as many cards as needed to specify the source values at the time 
points, allowing for 8 source values/data card, each being allotted 
10 columns of space (E) 

c) as many cards as needed to specify the time points, allowing for 


8 time points/data card, each being allotted 10 columns of space (E). 
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QUTPUT specifications: 
Only one type of output is allowed for any one network analysis, 


i.e., only one of NGRAPH, NALLOUT and JOUT can be nonzero. 
The following data cards follow immediately after the network description 
data cards. 
~ If NGRAPH # 0, we need NGRAPH (s 5) cards with the following information 
col 1-5, branch number (1) 
col 10, output type desired (I) 
1 = current desired 
O => voltage desired 
- If JOuT # 0, we need JOUT (= 200) cards with the following information 
col 1-5, branch number (I) 
col 10, output type desired (I) 
1 = current desired 
O = voltage desired 


- If NALLOUT # 0, no other data cards are required. 


Note that when one or more networks are analyzed in one batch, the 


last data card of the batch should be a blank card. 


oy) 


APPENDIX Cs 
” Sample Problems: 
Notation: 
. - circled numbers indicate node numbers 


uncireled numbers indicate branch numbers 


R's indicate resistances, in ohms 


- L's indicate inductances, in henrys 


t 


C's indicate capacitances, in farads 


NN is the number of nodes 


NB is the number of branches 


- H is the integration time step, in seconds 


TSTART is the starting time, in seconds 


TEND is the final time, in seconds 
- EPS desired performance criterion error 
Note, that only the pertinent computer output is included with the 


sample problems. 


Sample Problem #1: 
Butterworth Filter: 


ia 
I 


Q 
i] 





10 


= 0.01 


2x10 
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The input E is specified to be 


10 V 





time (seconds) 


0 2.5 x 10-6 5 x 10-6 


Specifications: 
- Initial time = 0 (TSTART = 0.0) 
- Final time = 1074 (TEND = 0.0001) 


{ime increment = 2 x 10° (H = 0.000002) 


- Desire error, at each time point, to fall below 10°’ (Eps 


tl 


107") 


1) 


- Desire results, at each time point, to be outputed (NSTEP 


- Desire tree and F to be outputed (NCONT = 1) 


Desire current to be in milliamps (SCALE = 10° ana NSCALE = 1) 


Desire the following graphical outputs 
Branch 1 voltage 
Branch 2 current (NGRAPH = 3) 


Branch 6 voltage 


Specify zero initial conditions. 
The required data cards, for this problem,may be seen on the following 
page. 


The central processor time for this problem was 2.954 seconds. 


DATA CARDS 
BBQUOUGDIEED 
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Sample Problem #2: 


Controlled Source DC problem: 





E =1 volt DC 

R, = RO = RE FLO 
NN = 5 

NB = 8 


Branch 3 is a current controlled current source, controlled by branch 4 
= x 
T, 5 I, 
Branch 8 is a voltage controlled current source, controlled by branch 7 
7 


Specifications: 
- since we have a resistive DC problem, we need the solution at only 
one time point 


set TSTART = 0.0 = TEND 


- desire error to fall below 107> 


(eps = 107) 


the time increment = 0 


(set H = 0.0) 


- desire results, at each time point, to be outputed (NSTEP = 1) 


do not desire tree and F to be outputted 


(NCONT = 0) 


desire current to be in amperes 


(SCALE = 1.0 and NSCALE = 1) 


desire all voltages and currents outputed 
(JOUT = 16) 
The required data cards, for this problem, may be seen on the 
following page. 


The central processor time, for this problem,was 0.24 seconds. 
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Sample Problem #3: 
D. C. amplifier: 





Transistor model: 


E 


}, 


The output will be the emitter current, and will be taken to be a 





current through a zero valued voltage source. 
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Modeling the controlled source in the way described in the section 


ft 


our network becomes 


on "Dependent Source Modeling,’ 


30 V 





Approximate theoretical gain: 
Ali5 


AT) 





Neglecting the 30K resistor, we can show that . 2.4, about a 


quiescent Lo of 1.1 ma. 


Specifications: 
- I, =A sin (1.8t) 
A = 0.1 milliamps 


t is in seconds 


initial time = 0 (TSTART = 0.0) 
final time =< 5 (TEND = 5.0) 


time increment = 0.1 (H = 0.1) 


Os dcetve tb 


proceed to next time point after 30 iterations (EPS = 107° ana 


desire error, at each time point, to drop to 10° 


NITE = 30) 


desire tree and F to be outputed (HCONT = 1) 


desire results, at each time point, to be outputed (NSTEP = 1) 
desire internal, automatic current scaling (NSCALE # 1) 


desire following graphical outputs 
Branch L current 
Branch & current 
(NGRAPH = 4) 
Branch 12 current 


Branch 6 current 


The required data cards for this problem may be seen on the following 


page. 
The central processor time, for this problem, was 8.24 seconds. 


DATA CARDS : PROBLEM #3 _S 
La] e] =f De] elephant of xp pclrbenpey Papp apap arrpepre 
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Sample Problem us: Compatible voltage source loop. 


=O 





R, = 22 


=1Q 
Rs 
NN = 4 
NB = 6 
Branch 3 is a current controlled voltage 


source, controlled by branch 6 


Vz = 3 * I 


Specifications: 


~ since we have a resistive DC problem, we need solution at only one 
time point. 


set TSTART = TEND = 0.0 
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- desire error to fall below 107? 


(EPS = 10°?) 


the time increment is arbitrary 


(set H = 0.0) 


desire results, at each time point, to be outputed 


(NSTEP = 1) 


desire tree and F to be outputed 
‘(NCONT = 1) 
- desire automatic internal scaling 


(NSCALE # 1) 


desire all branch currents and voltages to be outputed 


(JOUT = 16) 


The required data cards, for this problem, may be seen on the 


following page. 


The central processor time, for this problem, was 0.23 seconds. 
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Comments on Sample Problem #i: 

Note that when dependent and independent voltage source loops exist 
in a network, it is of paramount importance for them to be compatible, 
i.e., that Kirchhoff's voltage law be satisfied around that loop. 

Incompatible loops, of the above type, will yield incorrect executions 
of CANDO. The same holds for cutsets of independent and dependent current 


sources. 


TIME & 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
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BRANCK 
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BRANCH 
BRANCH 
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r wn fF &Y NN |= OBO WF F&F WwW N 
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CURRENT 
CURRENT 
CURRENT 
CURRENT 
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VOLTAGE 
VOLTAGE 


VOL TAGE: 


VOLTAGE 
VCLTAGE 


VOL TAGE. 


=,100006+01 
=,10000E+01 
oe 
~10000E401 
210000E¢0) 
210000E+01 
210900E+01 
.20000E40) 
63090001 
220000E+01 
210000E%01 
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Sample Problem #5: 
Four Step RC Ladder Network: 





R = 10° 
c =107? 
NN = 6 
NB = 9 


The input E is specified as follows 


t (seconds) 
0.0 0.01 
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Specifications: 
- initial time = O (TSTART = 0.0) 
- final time = 50 millisec (TEND = 0.05) 
~ time increment = 1 millisec (H = 0.001) 
- desire 10 iterative steps to be taken at each time point 
(EPS = 0.0, NITT = 10) 
- desire results, at each time point, to be outputed (NSTEP = 1) 
~ do not desire tree information and F to be outputed (NCONT = 0) 
- desire internal automatic current scaling (NSCALE # 1) 
- zero initial conditions 
- desire the following graphical outputs 
Branch 1 voltage 
} (NGRAPH = 2) 
Branch 5 voltage 
The central processor time, for this problem, was 8.26 seconds. 
The required data cards, for this problem, may be seen on the 


following page. 
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Sample Problem #6: 


Large time constant spread problem--quick response: 





We desire the step response associated with the small time constant. 


Specifications: 

- initial time = 0 (TSTART = 0.0) 

- final time = 50 microseconds (TEND = 5 x 10°”) 
:) 


- time increment = 1 microsecond (H = 10° 


- desire error to fall below 107) 


(eps = 10 2) 

- desire results, at each time point, to be outputed (NSTEP = 1) 
- do not desire tree information and F to be outputed (NCONT # 1) 
- desire internal, automatic current scaling (NSCALE # 1) 


- zero initial conditions 


- desire the following graphical outputs 
Branch 3 voltage 
Branch 4 current 
(NGRAPH = 1) 
Branch 2 current 


Branch 5 voltage 


The data cards, for this problem, may be seen on the following page. 


The central processor time, for this problem, was 4.11 seconds. 
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Comments on Sample Problem #6: 

The most important thing to notice in the computer output of this 
problem is that the voltage across Cs is, essentially, zero; i.e., the 
large capacitance is acting as a short circuit over the time interval of 
concern. 

The quick response steady state voltage across C, is 0.5 volts, and, as 
will be noted in problem #7, this value must be used as the initial condi- 


tion on C, when the slow time response is desired. 
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Sample Problem #7: 


Large time constant spread problem - slow response: 


The network is identical to that of problem #6. However, this time, 


we desire the step response associated with the large time constant. 


Specifications: 


initial time = 0 (TSTART = 0.0) 
final time = 500 sec (TEND = 500.0) 
time increment = 10 sec (H = 10.0) 


desire error to fall below 1074 or 20 iteration steps, per time 


point. (NITT = 20, EPs = 1072+) 

desire results, at each time point, to be outputed (NSTEP = 1) 
do not desire tree information and F to be outputed (NCONT = 0) 
desire internal, automatic current scaling (NSCALE # 1) 
initial voltages are 


V3 (9) =0.5V 


Vq5(0) =OV 


desire the following graphical outputs 
Branch 3 voltage 
Branch 4 current 
(NGRAPH = }) 
Branch 2 current 


Branch 5 voltage 


The data cards, for this problem, may be seen on the following page. 


The central processor time, for this problem, was 3.80 seconds. 
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Comments on Sample Problem #7: 

Note that the initial voltage on C3 has been set to 0.5 volts, which 
is the steady state value associated with the small time constant (see 
problem #6). Failing to do the above would yield incorrect results, since 
the trapezoidal integration scheme, applied to branch 3, would have a step 
size to time constant ratio of 10°, whereas a ratio of 0.1 is desirable 


for accurate integration. With the proper initial condition, C, becomes, 


3 
effectively, a voltage source, and the integration problem does not arise. 
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Sample Problem #8: 


Network with zero-valued elements: 


R, #1 NN = 5 
R, = 0 NB = 8 
R, = 2 
c, = 0 
Lg = 0 
CG, = 1 
Lg = 


Specifications: 


initial time =O (TSTART = 0.0) 

final time = 5 sec (TEND = 5.1) 

time increment = 0.1 sec (H = 0.1) 

desire tree and F to be outputed (NCONT = 1) 

desire error at each time point to fall below 107° (EPS = 107°) 
desire all voltages and currents outputed at every fiftieth time 


point (NSTEP = 50) 
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- desire internel automatic current scaling (NSCALE # 1) 

- zero initial conditions 

The data cards required for this problem may be seen on the following 
page. | 


The central processor time, for this problem, was 3.48 seconds. 











pero ToFOTOTooPSU SOT eToosToeovevarerorererersscorace cosa ocvoogOTaagsCOCCaae 


g,ANo> @# W3A180Ud . 


Comments on Sample Problem #8: 

Note that at the initial time point, all capacitances are tree branches 
and all inductances are links, as desired. For the remaining time points, 
zero valued R's and L's are tree branches and zero valued C's are links. 
Such 4 tree will always exist as long as there are no loops of short 
circuits (independent voltage sources, zero valued R's and L's), and no 
cutsets of open circuits (independent current sources and zero valued C's). 


result during the execution of the program. 


TImMe . 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
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BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 


BRANCH 
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06 
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06 


TIME = 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 
BRANCH 


- 
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VOI. TAGE 
VOLTAGE 
VO1, TAGE 
CURRENT 
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CURRENT 
CURRENT 
CURRENT 
CURRENT 
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010058E*02 
= 055976E°03 
= 056889E"01 
=.56889E"0) 
=210373E~02 
~0584R6E"01 
~.58486E"01 
7010058E*02 

°10058E*02 
-.29288E701 
= .28275E701 
= 613688E704 
~.28288E"01 
~.18666E*00 

010273E*%02 


APPENDIX D: 
CANDO FORTRAN IV LISTING 
(cpc 6400) 
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PROGRAM 


INTEGER 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
15) . 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
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CANDO (INPUT sOUTPUT) 


TYPE 
TEMP 
CONTYPEs 
SORTEMsC 
CNTSOR 
/BLOCK1/ 
/BLOCK2/ 
/BLOCK3/ 
/BLOCK4/ 
/BLOCK5/ 
/BLOCK6/ 
/BLOCK7/ 
/BLOCK8/ 
/BLOCK9/ 
/BLOCK10/ 
/BLOCK11/ 
/BLOCK12/ 


SORTYPE 
ONTTEM 


NNP( 200) »NP( 200) 
IBRAN( 200) sLEAV(200) sLENT( 200) 
TYPE (200) s VALUE (200) sF (499151) sNOEL (7) 
ITBRAN (200) sLEAVT( 200) sLENTT (200) 
IOUT (200) »NOUT (200) sITEST( 200) 
TEMP (7)sE( 200) sGRAD(200) »CONST (200) 
IDIMEN sAL PHA sFUNCT »NIC sNI TER sHs EPS »NNo NB 
V(20092)9C(20092) 
SORTYPE(200) sCONTYPE( 200) »KONBRAN(200) 
COND (200) sIDEL » TSTART » TEND »NSTEP 
SORVAL (59100) » TIMEPT (59100) sSNSOID (593) sNNI(5) 
CONTTEM(200) sSORTEM( 200) sCONDTEM( 200) sKONTEM( 200) 


/BK13/ ISTEP»OLDVAL (5) sSECVAL(5) sOLDTIME(5) sSECTIME(5) »NNU( 


/BLOCK14/ 
/BLOCK15/ 
/BLOCK16/ 
/BLOCK17/ 
/BLOCK18/ 
/BLOCK19/ 
/BLOCK20/ 
/BLOCK21/ 
/BLOCK23/ 
/BLOCK24/ 
/BLOCK25/ 


KP(5) 

NIT» JOUT »NGRAPH sNALLOUT 
GRAF (20095) sJGRAPH(5) s IGRAPH(5) sNPRINTs SCALE sNITT 
ITNsS(200) 

VAROUT (200) 

CNTSOR (200) 

NCONT »KLOOP »KPP 

NST(5) 
TEMPAR (200) sDUMMY (200) 
SCAL »NSCALE 
X(200)2¥(200) 


DIMENSION HH( 20920) 


NDATA=0 
NLNEAR= 
NIT=0 
ITN=1 
NPRINT= 
NIC=0 


0 


0 


IF (NDATAsEQe0) GO TO 2 


CALL SE 
TO=T=TU 
PRINT 1 
CALL SE 
CALL RE 
IDEL=NS 
NDATA=N 


COND (T) 


6s TO 
COND (TU) 
ADIN 
TEP=1 
DATA+1 


CALL KONST 


CALL IN 
GO TO 9 
NIC=NIC 


CREM 


+1 


CALL INCREM 
IF (NOEL(7)¢EQe0) GO TO 5 


NK=NB-N 
DO 4 I= 


OEL(7)+1 
NK »NB 


IF (SORTYPE(T)eEQe1HC) GO TO 4 


PPrPPrPTrrprrrrprrrprrrprrrrprrprrprrrrrrprrrrprprprrrrprprrrprprrrprprrrpypr > 


WOON OM FWNM & 


oe 


10 
ll 


12 


13 


14 


15 


16 


Cl 192)=C( 192) *SCALE 
CONTINUE 

NIT=0 

ITN=1 

IF (NLNEAReNE*O0) GO TO 6 

IF (NICeNEol) GO TO ll 

CALL PT 

IF (KPP.2EQ.0) GO TO 7 

GO TO 8 

SCAL=SCALE 

IF (NICeEQ@e1) CALL KONST 

GO TO ill 

CALL FCSM 

CALL KONST 

IF (NICeEQe0) GO TO 11 

DO 10 I=NNeNB 

C(Is2)=ClI92) *SCALE/SCAL 
C(191)=C(Il91)*SCALE/SCAL 
CONTINUE 

CALL CALCAL 

IF (NITTeEQe0) GO TO 12 

IF (NIT=NITT) 12914914 

CALL ERROR 

IF (FUNCT-EPS) 14914913 
CALL FNDGRAD 

CALL FLPOW (HHeNBsALPHA) 
ITN=ITN4+1 

NIT=NIT4#1 

CALL FNDALPH 

CALL CNGVARS 

GO TO ll 

IDEL=IDEL+1 

IF (IDELeNEeNSTEP) GO TO 15 
NPRINT=NPRINT+#+1 

IF (NGRAPHeGEel) CALL GRAPH 
IF (NALLOUTe¢GEol) CALL ALLOUT 
IF (JOUTeGEe1) CALL READOUT 
IDEL=0 

IF (NGRAPHeEQ.*201) GO TO 1 
IF (NIC=NITER) 39191 


FORMAT (1H1s20Xs*THE CENTRAL PROCESSOR TIME FOR THIS PROBLEM 
1E15e59* SECONDS*) 
END 
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SUBROUTINE READIN 
NB IS THE NUMBER OF BRANCHES 
NN IS THE NUMBER OF NODES 


INTEGER TYPE 

INTEGER TEMP 

INTEGER CONTYPE sSORTYPE 

INTEGER SORTEMsCONTTEM 

INTEGER TYTEMP 

COMMON /BLOCK1/ NNP( 200) »NP(200) 

COMMON /BLOCK2/ IBRAN(200) sLEAV(200) sLENT( 200) 

COMMON /BLOCK3/ TYPE( 200) sVALUE (200) 9F (499151) »NOEL(7) 

COMMON /BLOCK4/ ITBRAN(200) sLEAVT( 200) sLENTT (200) 

COMMON /BLOCK5/ IOUT(200)»NOUT( 200) sITEST(200) 

COMMON /BLOCK6/ TEMP(7)sITBB(200) »sTYTEMP( 200) sVALTEM( 200) 

COMMON /BLOCK7/ IDIMENsALPHA sFUNCT sNICoNITER sHs EPS sNNONB 

COMMON /BLOCK9/ SORTYPE(200) sCONTYPE( 200) »KONBRAN (200) 

COMMON /BLOCK10/ COND( 200) sIDEL »TSTART » TEND »sNSTEP 

COMMON /BLOCK11/ SORVAL (59100) s TIMEPT (59100) sSNSOID(593) sNNI(5) 
COMMON /BLOCK12/ CONTTEM(200) »SORTEM( 200) »CONDTEM( 200) sKONTEM( 200) 
COMMON /8K13/ ISTEPsOLDVAL (5) »SECVAL(5) sOLDTIME (5) sSECTIME(5) »NNUC 
15) 

COMMON /BLOCK15/ NIT»JOUT »NGRAPH sNALLOUT 


COMMON /BLOCK16/ GRAF (20095) sJGRAPH(5) » IGRAPH(5) sNPRINTsSCALE »NITT 


COMMON /BLOCK20/ NCONT sKLOOP 

COMMON /BLOCK24/ SCAL sNSCALE 

ISTEP=0 

IJK1=0 

I JK2=0 

KSOR=0 

KCURR=6 

READ 18 NNeNBsTSTART » TEND sHsNSTEP sEPS sNCONT 

IF (NNeEQeO) STOP 

READ 1959 NGRAPHsNALLOUT » JOUT sSCALE sNITT sNSCALE 
PRINT 209 EPSsTSTARTs TEND sH 

PRINT 219 NN»oNBsSCALE sNITT sNSCALE 

PRINT 22+ NSTEP 

IF (HeEQ@e00e0) GO TO 2 

NITER=(TENO-TSTART)/H 

GO TO 3 

NITER=0 

DO 9 T=1s5NB 

READ 239 TYPE(1) sIBRAN(1) »SORTYPE( 1) sCONTYPE(1) sKONBRAN(I) sLEAV(T) 
1lsLENT( I) sNCARDSsVALUE( 1) »sCOND(I) 

IF (NCARDSeEQ@e0) GO TO 9 

IF (TYPE(T)¢EQe1HJ) GO TO 4 

IF (SORTYPE( 1) eEQe1HSeOReSORTYPE(I) eEQe1HE) GO TO 5 
IF (SORTYPE(1)eEQ.1HPeORsSORTYPE(I)eEQe1HT) GO TO 6 
ISTEP=ISTEP+41 

IF (SORTYPE(CT) eEQ@e1HPsOReSORTYPE(I) ¢EQelHT} GO TO 8 
IF (SORTYPE(T) oEQe 1HSeOReSORTYPE(I) eEQe1HE) GO TO 7 


GO TO 9 
KSOR=KSOR+1 
NNJ(KSOR)=0 


READ 24 (SNSOID(KSOR»J) »J=193) 


FORDUDDDDOMDMAMDVMDDADVDDODDDODDODDVAVOTVOTDAOOVDDVDVDOTHOBDOWDWBODDOWOBBOBDWOOD 
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10 


ll 
12 


13 


14 


15 
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ow 


GO TO 9 

KSOR=KSOR+1 

NNJ(KSOR) =0 

READ 259 NNI(KSOR) 

NNN=NNI(KSOR) 

READ 269 (SORVAL(KSOR 9J) sJ=1sNNN) 
READ 269 (TIMEPT(KSOR»sJ) sJ=1 »NNN) 
GO TO 9 

KCURR=KCURR~-1 

NNJ(KCURR) =0 

READ 24 (SNSOID(KCURR »J) sJ=193) 
GO TO 9 

KCURR=KCURR=1 

NNJCKCURR )=0 

READ 259 NNI(KCURR) 

NNN=NNI(KCURR ) 

READ 26% (SORVAL(KCURR 9J) »J=1sNNN) 
READ 269 (TIMEPT(KCURRsJ) »J=1sNNN) 
CONTINUE 

IF (NGRAPHeEQe0) GO TO 10 


READ 27s (CIGRAPH(1) »s JGRAPH(1))»I=1eNGRAPH) 
IF (NALLOUT eGEeleOReNGRAPHeGE.1) GO TO 11 
READ 275 (CTOUT(I) sITEST(I)) sI=1.JOUT) 


PRINT 28 
PRINT 29 
DO 12 I=1sNB 


PRINT 309 IBRAN(I) sTYPE( 1) sVALUE(I) sLEAV(I) sLENT(I) sCOND(T) 


CONTINUE 

DO 14 I=1sNB 

IF (SORTYPE(1)¢NEe1H ) GO TO 13 
GO TO 14 

IF (TJK1eEQe0) PRINT 31 

TJK1=1 

A=7HVOL TAGE 

IF (TYPE(T) eEQelHJ) Aa THCURRENT 
PRINT 32» AP eras Poort cee 
CONTINUE 

DO 17 I319NB 

IF (CONTYPECI)eNEolH ) GO TO 15_ 
GO TO 17 

IF (1 JIK2eE&Q.0) PRINT 33 

TJK221 

IF (TYPE(T) sEQe1HI) GO TO 16 
A=7HCURRENT 


IF (CONTYPE(I)eEQe1HV) A=7HVOLTAGE 
. PRINT 349 AsIBRAN(1) sKONBRAN(T) 


GO TO 17 

A= 7HCURRENT 

IF (CONTYPE(T)eEQ@e1HV) A=7HVOLTAGE 
PRINT 359 AsIBRAN(T) sKONBRAN(I) 
CONTINUE 

IDIMEN=NB-NN+1 

KLOOP=0 

CALL PT 

IF (KLOOPeEQe1) GO TO 1 
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22 
23 
24 
25 
26 
27 
28 
29 
30 
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CALL FCSM 
RETURN 


FORMAT (159159E150e39E1503 9E15e3915sE15039I15) 

FORMAT (3159E15¢e5915915) 

FORMAT (1H19/////940X9* ERROR CRITERION = *9E15059///940X9% STARTI 
Sere es = *eE15 059% END TIME = *9E15¢59///s40X9*STEP SIZE = *» 

2 E1525) 

FORMAT (1HOs///s30Xs*NUMBER OF NODES = *sI59*% NUMBER OF BRANCHE 
1S = *9159///930X9*SCALE FACTOR = *9E615059///930X9* DESIRED NUMBER 
20F ITERATIONS AT EACH TIME POINT = *9159///930Xs*NSCALE = *915) 

FORMAT (1HOs///930X9*% NUMBER OF TIME ITERATIONS PER OUTPUT = *915) 

FORMAT (AlsI3sAlsAlsI391X9I291X91291X91292X%9E15049E1564) 

FORMAT ¢ 3E10¢3) 

FORMAT ( I5 } 

FORMAT ( 8E10¢3 ) 

FORMAT (215) 

FORMAT (1H1955Xs*THIS IS THE GIVEN NETWORK*) 

FORMAT (1HOs///955Xs*UNITS ARE OHMSs FARADS AND HENRYS *) 

FORMAT (1HOs 3X9*BRANCH NUMBER *92X91394Xe9* IS A * Al o* OF 
1 VALUE *s£12¢59* LEAVING NODE *sI30% AND ENTERING NODE *9139 
22X9* COND = *9E10¢3) 

FORMAT (1H19////930Xs*INDEPENDENT SOURCES*9///) 

FORMAT (1HOs10X9A7s * SOURCE BRANCH *9139* IS OF TYPE *sA1) 

FORMAT (1H=9///930Xs*CONTROLLED SOURCES *9///) 

FORMAT (1HOs10XsA7s* CONTROLLED VOLTAGE SOURCE#s1I39* IS CONTROLLED 
1 BY BRANCH*sI5) 

FORMAT (1HOs10X2A79* CONTROLLED CURRENT SOURCE*s1I39* IS CONTROLLED 
1 BY BRANCH*,15) 

END 


111 
1l2 
113 
114 
115 
116 
117 
118 
119 
120 
121 
‘122 
123 
124 
125 
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SUBROUTINE PT 

INTEGER TEMPO1 » TEMPO2 » TEMPOS sTEMPO4» TEMPOS sTEMPOT s TEMPO? 
INTEGER TYTEMP 

INTEGER TYPE 

INTEGER TEMP 

INTEGER CONTYPEsSORTYPE 

INTEGER SORTEMsCONTTEM 

INTEGER CNTSOR 

COMMON /BLOCK1/ NNP(200) »NP(200) | 

COMMON /BLOCK2/ IBRAN(200) sLEAV(200) sLENT( 200) 

COMMON /BLOCK3/ TYPE( 200) » VALUE (200) oF (499151) sNOEL(7) 
COMMON /BLOCK4/ ITBRAN(200) sLEAVT( 200) sLENTT (200) 

COMMON /BLOCK5/ IOUT( 200) sNOUT(200) sITEST( 200) 

COMMON /BLOCK6/ TEMP(7)»ITBB(200) s TYTEMP( 200) sVALTEM( 200) 
COMMON /BLOCK7/ IDIMEN+ALPHA FUNCT »NIC #NI TER Ho EPS »NNONB 
COMMON /BLOCK8/ V(20092)2C(20022) 

COMMON /BLOCK9/ SORTYPE(200) sCONTYPE( 200) eKONBRAN (200) 
COMMON /BLOCK10/ COND(200) sIDEL»TSTART»TENDsNSTEP 
COMMON /BLOCK12/ CONTTEM( 200) »SORTEM( 200) sCONDTEM( 200) sKONTEM( 200) 
COMMON /BLOCK15/ NIT »JOUT »NGRAPH sNALLOUT 

COMMON /BLOCK16/ GRAF (20095) sJGRAPH(5) sIGRAPH(5) sNPRINT»SCALE»NITT 
COMMON /BLOCK18/ CONST(200) 

COMMON /BLOCK19/ ‘CNTSOR(200) 

COMMON /BLOCK20/ NCONT »KLOOP sKPP 

COMMON '/BLOCK21/ NST(5) 

KK=NN=1 

KM=0 

KKTT=0 

IF (NICeNEeO) GO TO 5 

TEMP(1)=1HE 

TEMP(2)=1HV 

TEMP (3)=1HC 

TEMP (4)=1HR 

TEMP(5)=2HL 

TEMP(6)=1HI 

TEMP(7)=1HJ 

DO 3 K=1e7 

KT=0 

DO 2 I=1sNB 

IF (TYPE(I) sEQeTEMP(K)) GO TO 1 

GO TO 2 

KM=KM+1 

ITBRAN(KM) =I 

TYTEMP (KM) =TEMP(K) 

VALTEM(KM) VALUE (1) 

KT=KT+1 

CONTTEM(KM) =CONTYPE(T) 


SORTEM(K 
CONDTEM( 
KONTEM(K 
CONTINUE 
NOEL(K}= 
CONTINUE 
DO 4 [=] 
K=ITBRAN 


M)sSORTYPE (1) 
KM) =COND(f) 
M) =KONBRAN (1) 


KT 


oNB 
(1) 
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LFAYT(T)=LFAV(K) 

LENTTCI)=SLENT (CK) 

TYPFC(TIRTYTEMP(T) 

VALUE CT)=VALTEM(T) 
CONTYPELT)=SCONTTFM(T) 

SORTYPE( T)=SORTEM(T) 

COND (CT) =CONDTEM(T) 
KONBRAN(T)=KONTFM(T) 

CONTINUF 

TF (NOEL (4) 0¢FQ0e%) GO TO 15 
N=NOFLO1)4+NOFL(02)4NOEL (3) 41 
MMEN+NOFL(4)-1 

GO TO 6 

CONTINUF 

MM=ENRB—NOEL (7) 

N=NOFL(1)41 

M=MM—] 

DO 7 T=N5MM 

IF (TYPE CT) eFQe1HV). CONST(1T)=1-0F-50 
TF (TYPFCT) eFQe1HO) 69 TO B 

TF (TYPE (CT) eFQe1HR) CONST(T)=VALUF(T) 
IF (TYPE (CT) ¢FQelHL) CONST( I) =2eO#VALUECI)/H 
IF (TYPF(T) eFQe1HI) CONST(1T)=1-20F+50 
CONTINUE 

G9 TO 9 

TF (VALUF CT) ¢FQe9.60) CONST(T)=1.¢0F+51 
IF (VALUE CE) eNF e099) CONST(I) FH/(2 O#VALUF(T)) 
GO TO 7 

KPP=0 

IF (NeGTeM) GO TO 15 

MO 14 I=NoM 

KP=0 

AMIN=CONST(T) 

K=1+1 

MO 11 J=KoMM 

IF (CONST( J) eLTeAMIN) GO TO IN 

GO TO ll 

AMIN=CONST (J) 


— NP(T)}SJ 


KP=] 

KPD=} 

CONTINUF 

IF (KPeFO09) GO TO 14 
J=NP(T) : 


“TEMPOL=ITBRAN(T) 


TEMPO2=LEAVT(1) 
TEMPO3=LENTT(T) 
TEMPO4S=TYPE(T) 
TEMPN5=VALUF(1) 
TFMPNG=CONTYPF(T) 
TEMPOT=SORTYPF(T) 
TEMPOB=COND(1) 
TEMPO9=KONBRAN(T) 
TEMP1LO=CONST(I} 

IF (NICe~Q00) GO TY 12 
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» 12 


13 


14 
15 


16 


17 


19 
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TEMPL1=V(1T 9?) 

TFEMP12=C (19?) 
TEMP13=V(I91) 

TEMP14=C (191) 
ITBRANCTISITRRAN( J) 
LFAVT(CI)=LEAVT (J) 
LENTTCI)SLENTT(J) 
TYPEC(I)=TYPE(J) 
VALUF CT) =VALITE (J) 
CONTYPE (TY =SCONTYPE (J) 
SORTYPE(T)=SORTYPE (UY 
COND (T)=CONN (J) 
KONBRAN(T) SKONBRAN(J) 
CONST(CT)=CONST(J) 

TE (NTCeEQ20) GO TO 13 
V(Ie2)=V(Je2) 
C(192)=ClJe2) 
V(Ts1l)2#VJe1) 
C(I91)=C(J91) 
ITSRAN(J)=TEMPO1 
LEAVT(J)=TEMPO2 
LENTT(J)=TEMPO3 
TYPE(J)=TEMPO4 
VALUF (J) =TEMPOS 
CONTYPE(J)=TEMPN6 
SORTYPE (J) =TEMPY7 

COND (J) =TFEMPO8B : 
KONBRAN(J)=TEMPOS . 
CONST(J)=TEMP10 

Tf (NICSFEQ00) GO TO 14 
V(J92)STEMP11 
C(Js2)=TEMP12 
Vt(Jd0el)=TEMP13 
Cl(Jsl)=TEMP14 
CONTINUE : 
IF (2¢GTeKK) GO TO 32 
NO 25 T=29KK 
NE=LENTTC(T) 
NL=LEAVT(T) 

NP(1)=I 

NNP(T)=1 

MaT<=] 

MT=0 

KT=1 

DO 19 JJJ=1—M 
J=M4+1—“JII : 
IF (MTeFQeJ) GO TO 18 
IF (LEAVT(J)eFQeNEZE) GO TO 23 
IF (LENTT(CU)e¢EGQeNE) GO TO 24 
CONTINUF , 
CONTINUE 

IF (KTeFQel) GO TO 25 
KA=NNP(KT) 

KR=NP (KT } 


M=KR—1 
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110 


(1ii 


112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 


. 126 


127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
144A 
145 
146 
147 
148 
149 
150 
151 
152 
153 


155 
156 
157 
158 
159 
160 
161 
162 
163 


7] 


29 


26 


a7 


22 
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TE (MATOQN) AO TO 91 
“TSND(KT-1) 

KT=KT=] 

JF (KASFQe-1) 9 T9 2? 
NO =LFAVT(K8) 

ay ™9 17 

KT=KT-1 

GQ TV 20 

NEELENTT (KRY 

GQ TT. 17 

TE (LENTT OU) eEGeNL) GV TO 76 
NOSLENTT (J) 

KTSKT+1 

NOCKT)=J 

NND(KT)}=] 

VTsJ 

A=zJ~] 

G9 T. 17 

TF CLEAVT (SD) eFQeNL) GN TO 26 
NFSt_FAVT(J) 

“T=KT +1 

NOCKET) SJ 

NNP (WT) =] 

“T=J 

M=T=] 

G2 TA 17 

COMTT MIE 

GO TO 32 

LAF=NOFL (1) 

IF (Tel Fel AF) GO TY 99 
TEMPYATETTARAN(T) 
TOMPNO2=LEFAVTC-T) 
TEMPOASLENTT OT) 
TEVPNG=TYPE(T) 
TEMPASZEVALIUIF (CT) 
TFMPOGECONTYPE (TJ) 
TEMPOTESORTYPF(T) 
TEYPNBECONN(T) 
TEMPOAQEYONRBRAN(T) 

TE ONTC eI) GO TD FH 
TEMBIT=EVET 9?) 
TOAD P=C (1?) 
TEMDT3=V(T9)) 
TFMPYL4G=EC(T 91) 

CONTINUE 

AL RENQm] 

TF (NICeNF LO) NLREMNR=NOSL (7)+1 
NO 30 JMN=TsMNLA 

JF = JMN+] 

TTREASC IMPSTTRRAMG( JPY 
LE AVTCINIFLTAVT( JF) 
LENTT( US) SLENTT (UE) 
TYPE (UN) =TYPF (JE) 
VALUE CUN)EVALUFE (UF) 
CONTYPFE (JN) SCONTYPE( JE) 


VAN ANAAAVAANNANANDNANANANDNAANNNANANAANNANDNNNNANANDANAANAAANANNANA 


164 
165 


17? 
ere? 
175 
176 
177 
178 
179 
189 
181 
182 
183 
184 
185 
186 
187 
183 
189 
199 
191 
192 
192 
194 
195 
196 
LOT 
198 
199 
200 
201 
29? 
203 
204 
205 
206 
207 
ace 
209 
210 
211 
21? 
272 
214 
I15 
214 
217 
218 


229 
an 


31 
32 


24 


34 
35 
36 


37 


38 


39 


‘40 
41 


SORTYPFE (UN) =SORTYPF( JF) 


“COND (JN) =COND (JE) 


KONBRAN( JN) =KONRRAN( JF) 
IF (NICeFQe0) GN TO 29 
V(INs2)=V(IJE92) 
C(INs 2) SC( JE 92) 
V(INegl)=V(JF el) 
C(JNs 1) =C (JF 91) 

CONTTNUF 

CONTINUF 

NZ=NR 

TF (NICeNFeO) NZENR-NOFL (7) 
ITRRAN(NZ)=TEMPO] 
LEAVT(NZ)=TFMPO? 

LENTT(ONZ )=TEY“PO3 

TYPE(NZ) =STEMPO4 
VALUE (NZ )= TEMPOS 

CON TYPE (NZ) =TEMPON6 
SORTYPE(NZ)=TEMPO7 
COND( NZ) =TEMPO8 
KONBRAN(NZ)=TEMPOS 

IF (NICeFQLO) GO TO 31 
V(NZ92)=TEMP11 
C(NZs2)=TEMP12 
VONZs 1) =TEMP13 
C(NZg1)=TFMO1G 

CONTINUE 

IF (KKTTeFQ21) GO TO 35 
GO TO 16 

CONTINUF 

NO 33 L=loKK © 

IF (TYPF(L JF Qe TEMP(7)) GO TO 53 
CONT INUF : 

TF (NTCeNF 00) GON TON 36 
DO 35 I=NNoeNB 
IF (TYPF(T) -EQeTEMP(7)) GO TO 34 
GQ TO 35 

KKTT=1 

GO TO 27 

CONTINUE 

DO 37 IT=lsNAR 

K=TTRRAN(T) 
NP(1)=IBRAN(K) 

CONTINUF 

DO 38 T=1lsNBR 

CNTSOR(1)=0 

CONTINUF 

IF (JOUTeFQe0) GO TO 41 
DO 40 I=1sJOUT 

K=0 

K=K+1 

IF (IOUT(I) eNEeNP(K)) GO TO 39 
NOUT( 1) =K 

CONTINUE : 

IF (NGRAPHsEQe 0) GO TO 44 
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219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 


237 
238 
239 
240 
241 


243 
244 
245 
246 


4? 


44 


45 
46 
A7 
43 
49 


59 


] 


53 
54 
55 
56 
57 


58 
59 


4° 
41 
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D9 43 T=] 9NGRAPY 
wed 

K=%+] 

TF (TGRAPH(T) eNTeNP(K)) GO TO 42 
NOT(T) SK 

CONTINUE 

MEMOFL(1)41 

MaNReNIEL (7) 

NA SN T=NgM 

TE (TYPF(T) ef Me1H1) 40 TO 45 
IF (TYPF(T) = QelHV) SA TO 45 
GO TO 50 

K=KONBRAN(T) 

L=9 

L=L+1 

TE (LeGTeNOFL(1)) GO TO 47 

IF (NP(L)e¢NFRK) GN TA 46 

69 TH 40 

L=NS-NOFL(7) 

L=L+} 

IF (NP(L)eNFeK) GO TO 48 
CNTSOR(T)=L 

CNTSOR(L)=I 

CONTINUF 

TF (NCONTSCE 943) GD TO 54 
PRINT 55 

PPINT 5659 ((ToNP(T))91=19KK) 
PRINT 575 (CT NPC) ) 9 T=NNgNR) 
PRINT 58 

YY 51 T=1.N8 

PRINT 59s NO(T)sTYPT(T) VALUE CT) sLFAVTOL) sLENTT(I1) s CONT (T) 
CONTINUE 

ROTYRN 

PPTNT 69 

KLAOP=1 

RETURN 

PRINT 61 

KLAOP=1 

RETURN 


FCRMAT (TH1L es 2° Xs #CIRPESPONDENCE QP TWEEN ORIGINAL TAPNLAGY AND rw 


TTOPOLOGY #9 ///) 


FORMAT (THO 93°Xs#TREF PRANCHH 514 95Xs*CORRFE SPONDS TO POARCH %y14) 
FORMAT (1HM933X9% LINK %52XoT495X9*CORREGPONDS TA APANCH %574) 
FORMAT (1H1) 

FORMAT (1HO5 3X9#8RANCH NUMBEP %92XoIT324Xo% IS A *9Alo% OF 


1 VALUF #6512.59% LFAVING NONT He TBe% AND CNTERING NONE eo T3e 
2?X9% COND = %5F19,3) 


FORMAT (LHOo39Xs# THERE TS A VOLTAGE LOOP TM ME TWORe sy 
FOOMAT (THN g2°XS¥ETHERE TO A CURRENT GOMROE CLITCET TN NETWOP YH) 
aaah) 


AVVANANANANVNVANANNNNDNNAANANNAANNNANDNDNANDNNNANNANNNANNANNNANAAANADAAN 


274 
275 
276 
277 
278 
279 
289 
781 
28? 
283 
284 
285 
286 
287 
288 
289 
290 
29) 
292 
293 
294 
295 
296 
297 
298 
299 
306 
301 
302 
303 
304 
305 
396 
307 
308 
305 
310 
311 
312 


314 


319 
320 
221 
3?2 
322, 
324 


SUBROUTINE FCSM 

INTEGER TYTEMP 

INTEGER TYPE 

INTEGER TEMP 

COMMON /BLOCK1/ NNP(200) sNP(200) 

COMMON /BLOCK3/ TYPE (200) »VALUE (200) 9F (499151) »NOEL(7) 
COMMON /BLOCK4/ ITBRAN(200) sLEAVT( 200) sLENTT (200) 
COMMON /BLOCK5/ IOUT(10)sITEST(10) sNOUT(10) 
COMMON /BLOCK6/ TEMP (7)s1TBB(200) »TYTEMP( 200) »VALTEM( 200) 
COMMON /BLOCK7/ IDIMENsAL:/PHA s FUNCT sNIC »NITER sHo EPS #NNoNB 
COMMON /BLOCK20/ NCONTsKLOOP 

SET MATRIX TO ZERO 

KK=NN-1 

NLB=NB-NN+1 

DO 1 I=1sNLB 

DO 1 J=1sKK 

F(Jst)=000 

CONTINUE 

DO 8 I=NNsNB 

NE=LENTT (1) 

NL=LEAVT (1) 

1O=I-NN+1 

NP STORES THE PREVIOUS TREE BRANCH 

NNP STORES ITS DIRECTION 

NP(1)51 

NNP(1)=1 

M=1 

KT=1 

DO 3 J=MsKK 

IF (JeEQeNP(KT)) GO TO 3 

IF (LEAVT(J) eEQeNE) GO TO 6 

IF (LENTT(J)eEQeNE) GO TO 7 

CONTINUE 

M=NP(KT) 41 

KA=NNP (KT) 

KB=NP(KT) 

F(KB9T0)=000 

KT=KT=1 

IF (MeGEeNN) GO TO 4 

IF (KAeEQe1) GO TO 5 

NE=LENTT (KB) 


GO TO 2 

NESLEAVT (KB) 

GO TO 2 

F{(JsIO)=-1.0 

KT=KT+1 


IF (LENTT(J) sEQe NL) GO TO 8 
NE=LENTT (J) 

NP(KT)=J 

NNP(KT)=1 

M=1 

GO TO 2° 

F(JsTO}=100 

KT=KT+1 

IF (LEAVT(J)eEQeNL) GO TO 8 


900000 vovuDTVDUVUCU OUND OVUVVVOVUNUNUVVNUUOGVUNVNVVNOCVNOVONUONUNIITGIIGIVI VS... 


NNNMNN NNDB PPR ee ee ee 
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29 


19 


11 


12 


13 


15 


17> 


18 
19 


C23 


NE=LEAVT (J) 
NP(KT)=J 


NNP (KT) 
M=1 
GO TO 2 


ace 


CONTINUE 
IF (NCONT*EQe0) GO TO 16 
NLB=NB-NN+1 : 


KK=NN-1 


NCOUNT=—-1 


NT=NLB- 
NPAGF=1 


PRINT 179 NPAGE 
PRINT 18 
PRINT “199 (JsJ=19NLB) 


NT=NT-2 


5 


NCOUNT=NCOUNT+1 


NS1=1+2 


5¥NCOUNT 


NS2=25+2.5*NCOUNT 

IF (NTeLEeQ) GO TO 11 

DO 10 I=19KK 

PRINT 209 Is(F(IsJ)sJ=NS19NS2) 
CONTINUE 
NPAGE=NPAGE+1 


GO TO 9 


IF (NTe¢EQe0) GO TO 12 
GO TO 13 


NS3=1°. 
NS4=25 


Go TO 14 
NT=NT+25 


NS3=NS1 


NS4=NSL4NT=1 
CONTINUE 


pO 15 I 


=1eKK 


PRINT 205 Is(F(I9J)sJ=NS3 9NS4) 
CONTINUE 


RETURN 


FORMAT 


(1H1s40X9# F PORTION OF FUNDAMENTAL CUTS 


LF%sI5s//) 


FORMAT 
FORMAT 
FORMAT 
END 


(///) ; 
(1H0s7Xs25(2X91291X)) 
(1HOs1393X%s25(2X9F360)) 


13 


ET MATRIX *96X9*PAG 
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SUBROUTINE KONST 
CALCULATION OF CONSTANTS RELATED TO FLEMENT VALUES FOR USE IN BRAN 
RELATIONS AND GRADIENT 

INTEGER TYPE 
INTEGER TEMP 
INTEGER CONTYPE sSORTYPE 


COMMON 
COMMON 
COON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
15) 


/BLOCK3/ 
/BLOCK6/ 
/BLOCKT/ 
/BLOCK B/ 
/BLOCKOS 
/BLOCK10/ 
/3LOCK1I/ 
/BK13/ IS 
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TYPE( 200) sVALUIE (200) 9F (495151) sNOEL(7) 

TEMP (7) E200) sGRAD(200) sCONST (200) 

LD THEM s ALPHA sFUNCT sNICoNITERSHSEPS 9NNONB 
V(29092)9C(29992) 

SORTYPF (209) sCONTYPE (299) sKONBPAN( 200) 

COND (200) sIDELs TSTART s TEND sNSTEP 

SORVAL (52100) sTIMEPT (59109) »SNSOIDN (593) »NNIC5) 
TEP sOLDVAL (5) sSECVAL(5) sOLDTIME (5) »SECTIME(5) »NNJ( 


COMMON /3LOCK14/ NPD) 


COMMON /BLOCK.16/ GRAF (20025) sJGRAPH(5) s IGRAPH(S) »sNPRINTsSCALE®NITT 


COMMON /RBLOCK24/ SCAL »NSCALE 

IF (NICeNE29) GO TY 2 

No 1 12195 ; 

SECTIME(I)=TIMEPT( 192) 
OLDOTIME(1)=TIMEPT( 191) 
SECVAL( 1) =SORVAL(I 92) 
OLDVAL(1)=SORVAL(T91) 

NP(1)=2 

CONTINUE 

SCAL=SCALE 

IF (NSCALEeFilel) GO TO 5 

Kis 

DC 4 IT=NNeNB 

IF (TYPE( I) e«EQe1HR) GO TO 3 

GO TO 4 

IF (KWelQe0) SCALFEVALUF (I) 

KW] 

IF (VALUECT) eLTeSCALF) SCALE=VALUE(T) 
CONTINUE 

IF (KWeEQel) GO TO 5, 

DO 135 T=lsNN 

IF (TYPE(T)¢EQe1HR) GO TO 136 

GO TO 135 

IF (KWeFQeQ) SCALFRVALUE (I) 

KWwe] 

IF (VALUE(1)«GTeSCALF) SCALE=VALUF (IT) 
CONTINUE . 

IF (KWeFQe0) SCALT=1.0 

IF (SCALE eEQe00e3) SCALE=100 

N=NN-1 

NO 12 I=lsN 

IF (TYPE(I)e&Qe1LH=) GO TO 
IF (TYPECT) eFQelHV) 69 TO 
IF (TYPF(T) ef Qe1HC) GD TO 
IF (TYPE( IT) eEQelHR) GO TO 
IF (TYPE(T)eFQelHL) GO TO 
IF (TYPE(1)e&Qel1HI) 50 TO 
GO TO 12 


em OO NO 


He oO 


AM AM MMM AMM MMMM MMMM MMMM TMM AmMMaM m- 


P>Pryp 
900 
900 


eConoaourrwn 


TR CNET C ef ted) GSO TY 12 
VOETs2)=CONN(T) 

AAT. 12 
CONST CL) =VALNE(T) 


PP (CONTYPE CT) e& Ge THT) CONST CL) SCONSTOCL)/SSCOALS 


TF (MNTONE 20) 50 TO 1? 

VC Pe?) SCONN(T) 

GO TO 2 

COMET CLP aS OP eORVALUE CTY FSCALL 
IF OCNTCeRM.O) SO TO 12 
VtTA2PECONIOTS 

Bl PE. 1? 

COUSTOCLVSVALVECT IS SCALE 

TR (NECA .O) SO TO 12 
V(Ts2)2COND(T) 

af TO |? 

CONST CLP ES? eQFVALTIO CT) SCRE SCALL } 
PRONTO OME.) GD TN 12 
V(OTs2}49.% 

a: TO 12 

CONSTCDYSVALUECLE) 


TF CCONTY UOT oF ite LHV) CORST CPP SCONST C1) #5 


TF (NIC.NE 9) 35 TO 12 
ViTS2)20s 

CONT TMYE 

Mr 1G Laity NR 

IF (TYPE CT) oF Oe LI) TA TA 13 
TF (TYPF OL eh Qeltiv) C2 TO 14 
TF (TYPE CLI EQ.1LHT) 99 TO 15 
EF (TYPT C1 oie LHR) IC To 16 
[IF (TYP CLIT Q.1HL) GO T9 17 
IF (TYPPCL) EF Qe1HI} S56 TO 1B 
ot T3 1 9 

TF COMIC OMP.9) SO TS 1? 

OU Te SET A ee a 

so 79°19 

COUSTOULISVALHE CL) 


IF (CONT YP ER TD ete ltl + COMET OL ISCONST CI) /TCALE 


IF (MTCeN7.0) GO TO 19 
C(Ts?P)=9e% 

ae Ty Jo 

CONST CT) &2 OF VALUE (CT) 60 ALE /H 
it (Me Mego) Gn. pa 44 

aa ee eran 

Gey 3 
COMSTCL)=1eO/VALUE CL) ¥SC ALP 
TF (NICeN' 20) GO TA 19 
CllsP=CONDOT)FSCALE 

GQ TI 19 

CONST CL) SHSCALES( Pe OFVALUIE CTY) 
IF (NICANT.9) 59 TA 19 

CTs 2) = CONIC 1) SCALE 

59 TO 19 

CONST CL) =VALUECT) 


IF (CONTYPE (1) eEQe1HV) CONSTCT)SCONST CI) #S5CALE 


ace ee es 


AMM Meteo 


ebereG Mey in 


ee ees rey as 


PPT iy iy 


PTE Tr Sanna 


aM 


TERM ram in on 


‘t 


ao 


mim im: 


ft 
438 
49 
5 9) 


b2 
44 
A 
3D 
56 
57 
58 
5 
60 
&L 
be 
ae 
of 
oO 
ae a 
$6 
o7 
So 
Os 
i0 
7i 
72 
fas 
it 
rae 
76 
ye 


Vy 
3oO 
34 
G2 
bs 
34 
65 
36 
o7 
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ue 
2) 
Y1 
v2 
24 
Fa 
3% 
D6 
UT 


+ 


+3 
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- 102 


IF (NICeNE*0) GO TO 19 E 
C(1s2)=COND(1)*SCALE —°103 
CONTINUE E 104 
IF (NICsEQ.0) GO TO 20 E 105 
RETURN : E 106 
DO 21 I=lsNB | E 107 
IF (TYPE(1) eEQ@elHL) CONST(1)=020 E. 108 
IF (TYPE(T) e€Qe1HC) CONST(1)=000 E109 
CONTINUE E 110 
RETURN E 1ll 
END E 112- 


SUBROUTINE CALCAL 

INTEGER TYPE 

INTEGER CNTSOR 

COMMON /BLOCK3/ TYPE(200) sVALUE( 200) oF (499151) »sNOEL(7) 
COMMON /BLOCK8/ V(20092)9C(20032) 

COMMON /BLOCK7/ IDIMENsALPHAsFUNCT eNIC »NI TER sHe EPS sNNoNB 
COMMON /BLOCK19/ CNTSOR(200) 

N=NN=-1 

DO 1 I=loeN 

C(192)=000 

DO 1 J=NNeNB 

CCl s2)=Clls2)—-FUL sJ-N) #C (92) 

CONTINUE 

DO 2 I=NNeNB 

VtéIs2)=020 

DO 2 J=lsN 

VETs2.HVI 192) 4+F (Se I-N) #VEJ92) 

CONTINUE 

RETURN 

END 
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SUBROUTINE INCREM 

INTEGER TYPE 

INTEGER TEMP 

INTEGER SORTYPE sCONTYPE 

COMMON /BLOCK3/ TYPE(200) »VALUE (200) sF (499151) sNOEL(7) 

COMMON /BLOCK6/ TEMP(7)2E( 200) »GRAD(200) sCONST (200) 

COMMON /BLOCK7/ IDIMENs ALPHA FUNCT sNICsNITERsHeEPS »NNoNB 

COMMON /BLOCK8/ V{20092) sC(20092) 

COMMON /BLOCK9/ SORTYPE(200) sCONTYPE( 200) sKONBRAN (200) 

COMMON /BLOCK10/ COND(200) sIDEL»TSTART » TEND sNSTEP 

COMMON /BLOCK11/ SORVAL (59100) sTIMEPT (59100) sSNSOID(593) »NNIC5) 
COMMON /BK13/ ISTEP»OLDVAL(5) »SECVAL(5) sOLDTIME(5) »SECTIME(5) »NNJ( 


15) 


COMMON /BLOCK14/ NP(5) 

LOGICAL BOOL 

RBOOL=eTRUEe 

DO 1 T=#15NB 

V¢ Tol) =V(192) 

CONTINUE 

DO 2 IT=15NB 

C(Ts1l)=ClI92) 

CONTINUE 

TIME=TSTART+NIC#H 
M=NOEL(1)+NOEL (7) 

K=0 

DO 14 I=l»M 

J=] 

IF (1eGEeNOEL(1)4+1) J=NB—-I+NOEL(1)+1 
IF (JeGTeleANDeBOOL) GO TO 3 

GO TO 4 

BOOL=eFALSEe 

K=S5-I STEP 

IF (SORTYPE(J)eEGe1HC) GO TO 14 
K=K+1 

IF (SORTYPE(J) eEGelHEeOReSORTYPE(J) eEQ@e1lHS) GO TO 7 
IF (SORTYPE(J)eEQe1HP) GO TO 11 
IF (TIMEeLEeSECTIME(K)) GO TO 5 
IF (NNI (CK) eLEeNP(K)) GO TO 12 
NP(K)=NP(K)41 

NPR=NP(K) 

OLDTIME(K)=SECTIME(K) 
SECTIME(K)=TIMEPT(KsNPR) 
OLDVAL(K) =SECVAL(K) 

SECVAL(K) =SORVAL (K »NPR ) 

GO TO 15 

IF (TeGEe NOEL(1)+1) GO TO 6 

Wed cOLEGE VALI) Lt SECUAL | C1 SOCEVAL (C1127 IMEZOLDTIMETR D7 ESECTINER 
1K)-OLDTIME(K) ) 

GO TO 14 

ClJ92) SOLDVAL (K) +(SECVAL (K)-OLDVAL (K) )*(TIME-OLDTIME(K))/( SECT IME ( 
1K)-OLDTIME(K) ) 

GO TO 14 

W=SNSOID(K92) 

T=SNSOID(K 93) 

IF (SORTYPE(J)e¢EQe1HS) GO TO 9 


, 
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TF (JeGTel) GO TO 8 

VC Je2)=SNSOID (Ke 1) #EXP (We TIME-T) 
GO TO 14 
C(Is2)=SNSOID (Ks 1) EXP (VE TINE-T) 


°GO TO 14 


TF (JeGTel) GO TO 10 

VC J=s2)=SNSOID (Kel) #SIN(WHTIME-T) 
GO TO 14 
ClCJIs2)=SNSOID (Ke lL) *#¥SINCWHTIME-T) 
Ge TO 14 

CALL PERIOD (KoeJ)} 

GO TO 14 

L=NNI(K) 

1F (JeGTel) GO TO 13 
V(d92)=SORVAL (KsL} 

69 TO 14 . 

C(J92) =SO0ORVAL (Koel) 

CONTINUE 

RETURN 

END 


9 


ADADADAAADDAADAADHAAAAAD - 


te 


SUBROUTINE PERIOD (KoJ) 

INTEGER TYPE 

COMMON /BLOCK3/ TYPE (200) » VALUE (200) 9F (499151) sNOEL(7) 

COMMON /BLOCK7/ ID IMENs ALPHA » FUNCT sNICsNITER sHs EPS 9 NNoNB 

COMMON /BLOCK8/ V(20092)9C(20092) 

COMMON /BLOCK10/ COND (200) sIDELsTSTART » TEND sNSTEP 

COMMON /BLOCK11/ SORVAL( 59100) s TIMEPT (59100) »SNSOID(593) »NNI(9) 
COMMON /BK13/ ISTEP» OLDVAL(5) sSECVAL(5) sOLDTIME (5) sSECTIME(5) »NNJ( 


15) 
COMMON /BLOCK14/ NP(5) 
NNN=NNI(K) 


TIME=TSTART+NIC#H 

T=TIME-NNJ(K) #(TIMEPT(K»NNN)-TSTART) 

IF (TeLEeSECTIME(K)) GO TO 2 

NP(K)=NP(K)41 

IF (NP(K)eGTeNNI(K)) GO TO 4 

NPR=NP(K) 

OLDTIME(K)=SSECTIME(K) 

SECTIME(K)=TIMEPT(KsNPR) 

OLDVAL(K) =SECVAL(K) 

SECVAL(K) =SORVAL (K » NPR) 

GO TO 6 

IF (JeGEeNOEL(1)+1) GO TO 3 

V (392) ZOLDVAL (K) +(SECVAL (K) -OLDVAL (K) )#(T=OLDTIME (K))/(SECTIME(K)~ 
LOLDTIME(K)) A, 
GO TO 5 
C(Js2) =OLDVAL (K) +(SECVAL (K) =OLDVAL (K) )#(T=OLDTIME (K)) / (SECT IME (K)~ 
LOLDTIME(K)) -_ 

GO TO 5 

NNJCK)=NNJ(K) +2 

SECTIME(K)STIMEPT(Ks2) 

OLDTIME(K)=TIMEPT (Kol) 

SECVAL(K) =SORVAL(K 92) 

OLDVAL (K)=SORVAL (K 91) 

NP(K)=2 

GO TO l 

RETURN 

END 
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SUBROUTINE ERROR 
INTEGER TYPE 


INTEGER TEMP 
INTEGER SORTYPEsCONTYPE 
INTEGER CNTSOR 


COMMON 
COMMON 


COMMON /BLOCK7/ IDIMEN»sALPHA sFUNCT oNICsNITER sHe EPS oNNoNB 


COMMON 
COMMON 
COMMON 


/BLOCK3/ TYPE( 200) »VALUE( 200) 9F (499151) sNOEL(7) 
/BLOCK6/ TEMP(7)»E(200) sGRAD(200) sCONST( 200) 


/BLOCK8/ V(200s2) sC(200s2) 
/BLOCK9/ SORTYPE(200) sCONTYPE( 200) sKONBRAN (200) 
/BLOCK19/ CNTSOR(200) 


M=NOEL(1)+1_ 

N=NB=NOEL (7) 

DO 13 I=MsN 

IF (TYPE(IT) ¢EQe1HV) GO TO 3 

IF (TYPE(1) e¢EQe1HC) GO TO 5 

IF (TYPE(IT) e¢EQe1lHR) GO TO 7 

IF (TYPE(1)¢EQe1HL) GO TO 8 

IF (TYPE(I)-EQe1HI) GO TO 1 
E(1)=00 

GO TO 13 

K=CNTSOR(T) 

IF (CONTYPE(I)eNEelHv) GO TO 2 
E(1)=—-C(192)4+V(K92)*CONSTCT) 
GO TO 13 

EC T)=-C(192)4C(K 92) *CONST(T) 
GO TO 13 

K=CNTSOR(T) 

IF (CONTYPE(I)eNEelHV) GO TO 4 
E( 1 )s-V0l92)4V(K92)*CONST(T) 
GO TO 13 

ECT )e=-Vi192)4+C(K92)*CONST(I) 
GO TO 13 

IF (NICeEQe0) GO TO 6 

IF (1eGEeNN) GO TO 10 

ECT ==—(V¢C1s2)-VC191))4(C0192)4+C(191))} *CONSTET) 
GO TO 13 

E(1)=0e0 

GO TO 13 

IF (1eGEeNN) GO TO 11 
E(1)=-V(192)+C(1s2)*CONST(T) 
GO TO 13 

IF (NICeEQ20) GO TO 9 

IF (JeGEeNN) GO TO 12 
ECD=—(V¢To2)4VCT 51) }+0C(192)-ClIs1))*CONST(T) 
GO TO 13 

E(1)=060 

GO TO 13 

EC T)2=—(CC1ls2)4+C( lol) )+(V0192)-V0Il91))*CONST(I) 
GO To 13 
Et(l)=-C(192)4+V0192)*CONST(1) 
GO TO 13 
E¢lpe=(CCle2)—-Cllo1))+(V0C192)4+V0191))*CONST(T) 
CONTINUE 

FUNCT=0e0 
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DO 14 I=MeN _ 
FUNCT=FUNCT+E (1) ##2 
CONTINUE ; 
FUNCT=FUNCT/2 
RETURN 

END 
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SUBROUTINE FNDGRAD 


INTEGER CNTSOR 

INTEGER TYPE 

INTEGER TEMP 

COMMON /BLOCK3/ TYPE (200) sVALUE (200) 9F (499151) »NOEL(7) 
COMMON /BLOCK6/ TEMP(7)sE(200) sGRAD(200) sCONST( 200) 
COMMON /BLOCK7/ IDIMENsALPHA sFUNCT sNICsNITER sHe EPS »NNoNB 
COMMON /BLOCK19/ CNTSOR(200) 

COMMON /BLOCK23/ TEMPAR(200) sDUMMY (200) 


10 


M=NOEL(1)41 © 
MM=NB=NOEL (7) 

N=NN-1 

DO 5 I=19N 

IF (TYPEC(I) eEQe1HE) GO TO 1 
IF (TYPE(I)eEQe1HI) GO TO 3 
IF (TYPE(I) e¢EQelHV) GO TO 4 
TEMPAR(I)=-E(1) 

DUMMY (I) =CONST(I)*E(T) 

GO To 5 

K=CNTSOR(I) 

IF (KeEQe0) GO TO 2: 
TEMPAR(1)=000 
DUMMY (1) =CONST(K)*E(K) 

GO TO 5 

TEMPAR(1)=000 

DUMMY (1)=000 

GO TO 5 

TEMPAR(1)=000 
DUMMY (I) =-E(1) 

GO To 5 

TEMPAR(I)=~E(1) 
DUMMY (1) =060 

CONTINUE 

DO 10 I=NNeNB 

IF (TYPE(T)eEQelHJ) GO TO 6 
IF (TYPE(1)eEQelHV) GO TO 8 
IF (TYPE(I)eEQe1HI) GO TO 9 
TEMPAR(I)=CONST(I)*E(T) 
DUMMY (I) =-EC1) 

GO TO 10 

K=CNTSOR(I) 

IF (KeEQe0) GO TO 7 
TEMPAR( 1) =CONST(K) #E(K) 
DUMMY (1)=060 

GO TO 10 

TEMPAR(.1)=020 

DUMMY (1) =000 

GO TO 10 

TEMPAR(I)=-E(1) 
DUMMY (1) =000 

GO TO 10 

TEMPAR(1)=000 
DUMMY (1) =-E(T) 

CONTINUE 

IF (NICeNEeO) GO TO 13 
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DO 12 I=MsMM 


IF (TYPE( 1) eEGelHCeOReTYPE( I} eEQel HL) GO TO ll 


GO TO 12 

TEMPAR(1)=009 

DUMMY (1) =020 

CONTINUE 

IF (MeGTeN) GO TO 18 

DO 14 I=MsN 

GRAD(1)=TEMPAR(I) | 

DO 14 J=NNeNB = 
GRAD(1)=GRAD(I)4+F (I 9J-N) #TEMPAR (J) 
IF (NNeGTeMM) GO TO 19 

DO 15 I=NNsMM 

GRAD(1)=DUMMY (1) 

DO 15 JeleN 
GRAD(I)=GRAD(1)-F(JsT-N) *DUMMY (J) 
IF (NICeEQe0) GO TO 16 

RETURN 

DO 17 I=lsNB 


IF (TYPE(1T) eEGelHC) GRAN(1)=000 
LF (TYPECIT) eEQelHL) GRAD(1)=000 
CONTINUE 

RETURN 

END 
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SUBROUTINE FNOALPH 
THIS SUBROUTINE FINDS THE DISTANCE WF SHOULD GO ALOND THE GRADIENT 
INTEGER TYPE 

INTEGER TEMP 

COMMON /BLOCK3/ TYPE(200) sVALUF (200) sF (4959151) sNOEL(7) 
COMMON /BLOCK6/ TEMP(7)sE( 200) sGRAD(209) sCONST( 200) 
COMMON /BLOCKT/ TDIMFENsALPHAsFUNCT sNIC oNITER sHe EPS oNNONR 
COMMON /BLOCKB8/ V(20092)9C( 2092) 

COMMON /BLOCK12/ CONTTEM( 290) sSORTFM( 200) sCONRTEM( 209) sKONTEM( 200) 
COMMON /BLOCK15/ NIT 

COMMON /BLOCK17/ ITNsS(200) 

FZERO=FUNCT 

M=NOFL(1)+1 

N=NB=NOEL( 7) 

MM=NN-=1 

PRON=ENeD 

IF (MeGTeMM) GO TO 5 

DO 1 I=MsMM 

PROD=PRON+GRAD(1)*S(1) 

CONTTEM( I) =V¢192) 

VCIs2)=V(192)-S(1) 

CONTINUE 

IF (NNeGTeN) GO TO 6 

NO 2 I=NNoN 

PROD=PROD+GRAD(1)*S(1) 

CONTTEM(1)=C( 192) 

C(192)=C(192)-S(1) 

CONTINUE 

CALL CALCAL 

CALL ERROR 

FGRAN=FUNCT 

IF (MeGTeMM) GO TO 7 

DO 3 1=M yMM 

VOIs2)=CONTTEM(1) 

CONTINUE 

IF (NNeGTeN) GO TO 8 

DO 4 IT=NNoN 

C(192)=CONTTEM(I) 

CONTINUF 

Al_PHA=PROD/( FGRAD+PRON=-FZERO) 

ALPHA=ALPHA/2 0 

RETURN 

END 
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SUBROUTINE CNGVARS 


THIS SUBROUTINE CHANGES THE CURRENTS ALONG THE GRADIENT BY THE 


CONSTANT ALPHA 

INTEGER TYPE 

INTEGER TEMP. 

COMMON /BLOCK3/ TYPE (200) sVALUE( 200) oF (499151) »NOEL(7) 
COMMON /BLOCK6/ TEMP (7) 9E(200) »GRAD(200) sCONST (200) 
COMMON /BLOCK7/ IDIMEN ALPHA » FUNCT sNIC »NITER sHo EPS oNNONB 
COMMON /BLOCK8/ V(20052)9C(2002) 

COMMON /BLOCK17/ ITNsS(200) 

M=NOEL(1)+1 

N=NB=NOEL(7) 

IEND=NN-1 

IF (MeGTelEND) GO TO 3 

DO 1 I=MsTEND , 

V( 192) =V¢192)-ALPHA*S(1) 

CONTINUE 

IF (NNeGTeN) RETURN 

DO 2 I=NNeN 

C(1s2)2C¢1 92) ~ALPHA*S(1) 

CONTINUE 

RETURN 

END 
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SUBROUTINE FLPOW (HH»NKsALPHA) 

INTEGER TEMP 

INTEGER TYPE 

COMMON /BLOCK3/ TYPE( 200) »VALUE (200) 9F (499151) »sNOEL(7) 
COMMON /BLOCK6/ TEMP (7)2E(200) sGRAD(200) sCONST(200) 
CUMMION /BLOCK17/ ITNsS(200) 

COMFION /BLOCK25/ GRAD2( 200) sDUMMY (200) 
DIMENSTON HH( NBs NR) 

IF. (ITNeGTel) GO TO 3. 

M=NOEL(1)41 

N=NB=-NOEL (7) 

DO 2 T=MeN 

00 1 JeMoN 

HH(I9J)=060 

HH(IsIT)=1-0 

GRADB( 1) =GRAD(1) 

S(I)=GRAD(1) 

RETURN 

GHG=06« 

SG=06. 

DO 4 IT=MsN 

SCI) =ALPHA*S(1) 

GRADB(1)=S RAD(1)-GRADA(1) 
SG=SG+S(1)*GRADB(T) 

DO 5 IT=MaN 

DUMMY(T)=06 

DO 5 J=MoN 
GHG=GHG+GRADB (1) *GRADB (IJ) *HHCI »J) 
DUMMY.CT) sDUMMY (1)+HH( TsJ)*GRANDE (J) 

DO 7 I=M>5N 

GRADB(1T)=0e0 

DO 6 J=HIeN 
HH(IsJ)=HHCIsJ)-(501) 4504) /SG)— (DUMMY (1) #DUMMY (J) /GHG) 
HH(JsI)=HH( IJ) 


DO 7 J=MsN 


GRADB(1)=GRADBC1T)+HH¢1»sJ)*GRAD(J) 
DO 8 I=MsN 

S(1)=GRADB(T) 

GRADBII)=GRAD (1) 

RETURN 

END 
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SUBROUTINE READOUT 

INTEGER TYPE 

INTEGER TEMP 

COMMON /BLOCK3/ TYPE(200) sVALUE (200) sF (499151) sNOEL(7) 
COMMON /BLOCK5/ IOUT(200) sNOUT( 200) sITEST(200) 

COMMON /BLOCK6/ TEMP(7)+E(200) sGRAD( 200) »CONST (200) 
COMMON /BLOCK7/ IDIMEN»sALPHA sFUNCT sNIC sNITER sHo EPS sNNoNB 
COMMON /BLOCK8/ V(20092)2C(20092) 

COMMON /BLOCK9/ SORTYPE( 200) sCONTYPE( 200) sKONBRAN (200) 


wa 
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COMMON 


/BLOCK10/ COND(200) sIDEL»sTSTART s TEND sNSTEP 


COMMON /BLOCK15/ NITsJOUT »NGRAPH sNALLOUT 


COMMON /BLOCK16/ GRAF (20095) sJGRAPH(5) » IGRAPH(5) »NPRINT»SCALE sNITT 


COMMON /BLOCK18/ VAROUT(200) 
N=NN=1 | 

DO 4 I=1sJOUT 

K=NOUT (1) 

IF (NOUT(I)eLTeNN) GO TO 2 
IF (ITEST(1)e&Qe1) GO TO 1 
VAROUT(1)=V(K92) 

GO TO 4 

VAROUT (1) =C(K 92) /SCALE 

GO TO 4 

IF (ITEST(1)eEQ00) GO TO 3 
VAROUT (1) =C(K 92) /SCALE 

GO TO 4 

VAROUT(1)=V(K92) 

CONTINUE 

TIMESH*NIC4+TSTART 

IF (NICeEQe0) PRINT 6 

PRINT 7s TIME 

DO 5 1=1sJOUT 

A=7HVOLTAGE 

IF (ITEST(1) eEQ@el) A=7HCURRENT 
PRINT 8» IOUT(I) sAsVAROUT (1) 
CONTINUE 

RETURN 


FORMAT (1H1) 
FORMAT (/////s10Xe*TIME = *9E1505) 


~~ FORMAT (1HOs10X»*BRANCH* 9 1595X9A7910X%9E1505) 


END 
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SUBROUTINE GkAPd 


INTESER 


ANP 


INTEGER TYPE. 


Cun SO 
COV AON 
TG AON 
COSAON 
COMA ON 
CPiAON 
COMMON 
COnriON 
Coa doh 





No 3 TF 


/BLOCAL/ ANPC2900) shP C200) 

(LOCK 3/ TYPE (290) eVALiFE 0290) 9F (499151) sNOFEL C7) 

/BLOCKG/ TTRRAN( 200) eLEAVT (200) sLENTT (C290) 

/BLOCKT/ IDTME NS ALPHA S FUNCT ANIC ONITER sHe EPS oNN ONE 

/BLOCKB/ V>20092?)9C(20092) 

/LOCKLO/ COND( 206) sIDEL» TSTART es TEND otISTEP 

/RLOCK15/ NIT JOUT sNGRAPH sNALLOUT 

/BLOCK16/ GRAF (20025) »JGRAPH(5) » IGRAPH(5) »NPRINTsSCALEeNITT 
/BLOCK21/ NST(5) 

1s NGRAPH 


K=ENST(T) 


IF (JGRAPN]) 
GRAF CANPRINT® 1D) 
Gd TO 3 


het) GD TO 2 
HBC (K 92) /E CALE 


GRAF (CNPRENTs TL) EVC% 92) 
CONT TAGE 
E (NPRINTe LT oe 2G AND? (NI TERANIC) GF eNSTEP) RETURN 


i 
yo 9 Ts 
PRINT 10 


1» NGRAPH 


A= 7HVOLT AGE 

LE (JGRAPH( 1) eEGel) ASTHCURRERT 
PRINT lis IGXAPHCT) 5A 
AMAXSGRAF (LT) 

AHAX#ASS CAvIAX ) 


Oo 4 J= 


2 NPISI ANT 


A=GRAF (Je) 

AZ=ABPS{(A) 

IF (AeGTeAMAX) AMAXBA 
CONTINJEF 


D0 5 J 


LoNPRINT 


NP (J) =GRAF (Jo 1)#59 eU/AHAX4+51 69 
CONTINUE 


DG 6 JF 


12101 


ANP( J) =]He 
CONT TNUF 
PRINT 12% (ANP(J)9J=19101) 


99 7 J 


12192 


ANP (J)=1H 

CONTINUE 

ANP (51) = Lite 

DO 8 €K=1L sNPRINT 
LL=NP(KK) 

IF (LL eGFe51) LL=LLt1 
ANP(LL)=1n+ 
TIMFATSTARTER STE PEER EHH 
PRINT 139 TIE sGRAF (KK 1) 9 CANP OS) 9 Jl 9102) 
ANP(LL) =1Lti 

CONTI NUF 

COMTINGE 


NGRAPHS 


RETURN 
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APPENDIX E 


General comments regarding use of CANDO. 


Ts 


Number of nodes cannot exceed 50 

Number of branches cannot exceed 200 

Cannot exceed 5 arbitrary or periodic independent sources, and 
cannot exceed 5 sinusoidal or exponential sources. There is no 
limit on the number of constant independent sources. We can have 
at most 100 time points to describe an arbitrary independent source. 
Cannot exceed 5 graphical outputs 

There is an internal limit of 200 output points per graph. 

Branch and node numbering is arbitrary, but branch numbers cannot 


exceed three digits and node numbers cannot exceed two digits. 


Core requirements 


Since using the Fletcher-Powell minimization algorithm requires the 


storing of an NB x NB matrix, it would, in general, be inconvenient to 


always dimension the above matrix 200 x 200. 


Thus, to save core, one need only replace the DIMENSION statement in 


the MAIN program by 


DIMENSION HH(N, N), X(N), X(N) 


where N is the maximum number of branches to be encountered in the network(s) 


to be solved. 


With the above dimensioning procedure, the core requirements for CANDO 


are, approximately, 


27000 + n° (decimal) 


storage locations. 


ww 
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There must be at least one branch which is not an independent 
source, in any given network. 

CANDO can handle zero valued R's, L's and C's as long as there 
are no loops of zero valued inductances or resistances (short 
circuits) and no cutsets of zero valued capacitances (open circuits). 
In this context, independent voltage sources are equivalent to 
short circuits and independent current sources are equivalent to 
open circuits. 

H, the time increment, cannot be zero for networks containing 
reactive elements (L's and C's). 

Zero is not a valid branch or node number. 

CANDO can be utilized for batch processing; simply stack data 
card sets, for the networks to be solved, one behind the other, 


ensuring that the last card in the stack is a blank card. 
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